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I  INTRODUCTION 


This  is  the  final  report  of  a  series  that  was  initiated  under  a 
former  contract  (AF  29(601  J-199-J ,  ARPA  Order  No.  6-58,  Task  5)  and  has 
continued  for  a  period  of  approximately  3-1  2  years.  This  report  contains 
much  of  tile  earlier  work,  brought  up  to  date,  as  well  as  new  material 
generated  during  the  final  phase  ol  the  present  contract. 

The  objective  of  this  study  has  been  to  learn  as  much  as  possible 
about  the  relaxation  of  the  lower  ionosphere  following  a  pulse  of 
ionizing  radiation. 

Ideally,  the  approach  would  be  to  measure  several  ionospheric 
parameters  directly,  simultaneously,  and  continuously  during  the  course 
of  an  ionospheric  disturbance.  The  parameters  of  interest  would  be 
radiation  flux,  electron  concentrations,  ion  species  and  concentrations, 
and  gas  and  electron  temperatures.  For  many  reasons,  which  include 
limitations  of  the  state  of  the  art  and  financial  support,  the  ideal 
experiments  are  not  likely  to  be  performed  in  the  near  future.  Therefore, 
the  approach  adopted  in  this  study  has  been  to  review  published  and 
inpubi i shed  measurements  of  the  above  parameters,  especially  those  made 
duting  disturbed  conditions.  Pertinent  geophysical  measurements  were 
considered,  together  with  laboratory  data  and  current  theory,  in  an 
effort  to  delineate  the  phenomena  involved  and  subsequently  to  estimate 
the  rates  of  the  various  relaxation  processes. 

The  report  is  divided  into  two  principal  parts — ionization  processes 
anil  de-ioni zation  processes.  This  separation,  as  the  reader  will  see, 
is  not  clear-cut.  Inasmuch  as  natural  ionization  processes  must,  at 
present,  be  studied  by  considering  fragmentary  observations  of  the 
ionizing  flux  together  with  the  resulting  radio  signal  absorption,  an 
inescapable  dependence  on  de-ionization  processes,  which  is  not  yet 
clearly  defined,  always  exists.  Hence,  in  Section  III  A,  a  model  of 
solar  flare  ionization  is  used  to  deduce  an  effective  dissociative 


recombination  coefficient;  the  other  necessary  relaxation  parameters 
were  obtained  from  the  literature.  In  III  B,  solar  X-ray  observations 
were  used  with  parameters  discussed  in  Section  IV  to  gain  an  insight 
of  quiescent  solar  effects.  In  III  C,  a  model  of  auroral  electron- 
bremsstrahl ung  events  is  constructed  utilizing  parameters  derived  in 
IV  to  check  the  adequacy  of  those  parameters  in  describing  ionospheric 
perturbations.  In  a  similar  vein,  the  associative  recombination  co¬ 
efficient  discussed  in  IV  C  was  derived  from  observations  of  solar 
flares,  and  sunrise-sunset  radio  effects  are  used  in  IV  A  to  evaluate 
possible  candidate  negative-ions.  Also  in  Section  IV  C,  the  dissociative 
recombination  coefficient  derived  from  solar  flare  studies  is  shown  to 
be  consistent  with  laboratory  and  other  ionospheric  measurements. 

Thus,  although  definitive  atmospheric  measurements  have  not  been 
made,  a  reasonable  consistency  was  found  to  occur  between  estimates, 
laboratory  measurements,  and  pertinent  geophysical  observations  and 
models.  However,  this  apparent  consistency  should  not  lull  the  reader 
into  complete  complacency.  The  inherent  errors  in  many  of  the  necessary 
assumptions  are  such  that  it  would  not  be  at  all  surprising  to  find  that 
a  good  measurement  of  one  or  more  of  the  important  parameters  will  show 
these  models  and  assumptions  to  be  far  from  adequate. 

Significant  data  have  been  and  are  being  acquired  that  have 
materially  increased  our  understanding  of  ionospheric  processes.  Present 
knowledge  is  far  from  complete,  of  course.  Although  significant  first 
steps  have  been  taken,  it  is  clear  that  the  entire  field  of  solar- 
terrestrial  relationships  is  only  beginning  to  develop.  A  continuing 
review  of  laboratory  and  geophysical  observations  is  certainly  required; 
definitive  in  situ  measurements  are  urgently  needed. 

Earlier  work  on  this  study  is  described  in  the  following  reports 
and  papers . 

i 1  Poppoff,  I.  G  ,  K.  C  Whitten,  and  R .  L  Ludwig.  The 
Determination  of  Ionospheric  Recombination  Coefficients 
Phase  I:  Feasibility  of  Utilizing  Solar  Hare-Sudden 
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Ionospheric  Disturbance  Relationships.  Final  Report 
Contract  No.  AF  29 (601 )-1994 ,  ARPA  Order  No.  6-58, 

Task  5.  GRD-TR-60-289 .  Stanford  Research  Institute, 

May  26,  1960. 

2  Poppoff,  I  C.  and  R.  C.  Whitten,  Determination  of  Ionos¬ 
pheric  Recombination  Coefficients,  Semiannual  Report  No.  1 
Contract  AF  la  v  „0-l  j-8355  .  Stanford  Research  Institute, 
September  30,  1961 

3  Whitten,  R  C  and  I  G  Poppoff,  A  Model  of  Solar  Flare- 

Induced  Ionization  in  the  D-region.  J.  Geophys.  Res.  66, 
2779  (1961;.  '  — ’ 

X-1  Whitten,  R.  C.  and  I  G  Poppoff  Associative  Detachment 

in  the  D- region.  J .  Geophys.  Res.  67,  1183  1962  . 

(o  j  Poppoff,  I.  G  and  R  C,  Vihitien.  D- region  Ionization  by 

Solar  X-rays.  J  Geophys.  Res.  67,  2986  1962} 

i®)  Poppoff,  I.  G  and  R  C  Whitten.  Determination  of  Ionos¬ 

pheric  Recombination  Coefficents.  Semiannual  Reports  2  and 
3.  Contract  AF  19  601  -8355.  Stanford  Research  Institute, 
September  30,  1962 
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SUMMARY  AND  CONCLUSIONS 

pertaining  to  ionization  and  dc-ionization 
in  the  lower  ionosphere  was  compiled  from  various  sources.  Part  of 
the  data  was  derived  by-Uie  authors  from  studies  of  geophysical  ob¬ 
servations,  part  of  it  from  published  and  unpublished  laboratory 
measurements,  and  part  from  other  ionospheric  studies.  Models  of  lower 
ionospheric  processes  were  prepared  which  are  consistent  with  presently 
available  data  from  these  sources 

\ 

The  following  conclusions  may' be  drawn  from  this  study: 

\ 

A.  Ionization  Processes 


II 


A 


A  body  of  information 


Solar  Flare-Induced  Ionization 


Solar  flare  effects  can  be  considered  to  be  X-ray 
ionization  phenomena.  Consideration  of  sudden  cosmic 
noise  absorption  events  leads  to  an  effective  dissocia¬ 
tive  recombination  coefficient  of  about  2  x  10-7  cm3  sec*1 


Ionization  by  Background  Solar  X-rays 

Normal  D-rcgion  ionization  sources  have  not  been  definitely 
established.  Consideration  of  present  knowledge  of  nitric 
oxide  recombination  coefficients  and  background  solar  X-ray 
fluxes  indicates  that  the  role  of  solar  X-rays  may  be  greater 
than  previously  supposed  and  that  the  role  of  solar  hydrogen 
Lyman— (X  radiations  may  be  less  important. 

Ionization  by  Auroral  Electrons 


Studies  of  auroral  e lec t ron-bremsst rah  lung 
that  at  least  the  temporal  characteristics 
radio  noise  absorption  are  consistent  with 
suggested  in  this  study.  The  model  is  not 
enough  to  judge  its  adequacy  for  computing 
of  the  absorption. 


events  indicate 
of  associated 
the  model 
yet  complete 
the  magnitude 
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B 


De-ionization  Processes 


1.  Negative  Ions 


a.  Significant  concentrations  of  0~  and  H~  ions  in 
the  lower  ionosphere  are  not  consistent  with 
observed  "sunrise-sunset"  effect.  0“  may  be 
important  in  the  upper  D-region 

b.  02  ,  OH  ,  03  ,  and  or  N02  may  be  the  dominant 
negative  ion  species  without  contradicting  present 
knowledge . 

c.  Radiative  attachment  to  02  and  three-body  attachment 
to  0  are  probably  not  important. 

d  Associative  detachment  is  probably  the  dominant 
nighttime  mode  of  electron  detachment  from  the 
species  02~  and  Oil-.  Detachment  from  N02-  and  Oa- 
can  occur  via  dissociative  processes. 

e.  Previously  proposed  values  of  photodetachment  rate 
for  02-  and,  perhaps,  o-  are  not  reliable. 

f.  Diffusion  to  dust  is  an  unimportant  mode  of  negative 
ion  removal . 

Charge  Transfer  and  Ion-Atom  Interchange 


The  following  rate  constants  were  derived 


(1) 

0+  +  02 

•  0  +  02+, 

k 

2 

X 

10-1 1  cm3  sec-1 

(2) 

n2+  +  0 

. n2  +  o', 

k 

2 

X 

10-1'  cm3  sec-1 

(3) 

n2  +  02 

*  N  2  +  O2 

+ 

k 

- 

2 

x  10-,°  cm3  sec 

(*») 

o+  +  n2 

•  N0+  +  N , 

k 

■S. 

6 

X 

10” *  3  cm3  sec- 1 

b.  Rate  constants  for  1  and  (3)  above  are  in  good 
agreement  with  recent  laboratory  measurements,  whereas 
the  rate  constant  for  (4)  is  not.  The  discrepancy 
cannot  be  explained  at  this  time. 

c.  The  dominant  ions  in  the  D-region  under  disturbed 
conditions  must  be  02+  and,  or  N0+ .  The  model  pre¬ 
sented  in  this  report  cannot  be  extended  below  ~  75  km 
because  of  the  possible  importance  of  such  ions  as 

N 3  ,  N„+,  and  03+ . 


G 


3. 


Dissociative  Recombination 


a.  Probable  values  of  important  dissociative  recombina¬ 
tion  coefficients  in  the  D“  and  E*  regions  are: 

(*)  ^'^2  )  =  5  x  10  s  T“3  4  cm3  sec-1 

(2)  aD(°2+)  =  9  x  10-5  T- 1  cm3  sec-1 

3)  Qp(NO  )  =  2.6  x  10-4  T“ 3  2  cm3  Sec-1 

b  Most  important  coefficients  in  the  upper  D-region 
and  the  E-region  are  Cyo/)  and  a^NO*). 

c.  Identities  of  dominant  positive  ions  in  the  lower 
D-region  are  not  known. 

4 ■  Ion-Ion  Recombination 

a.  Two-body  reactions  are  important  in  the  D-region. 

b.  Three-body  reactions  are  not  important  above  ~  -15  km 
but  may  be  important  in  low-altitude  effects  of  nuclear 
explosions . 

5-  Diffusion  to  Dust 


This  is  not  likely  to  be  important  in  the  normal 
unperturbed  ionosphere  or  during  perturbing  events. 
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Ill  IONIZATION  PROCESSES 


Solar  FI are-I nduced  Ionization 
1.  Introduction 


It  has  long  been  known  that  so-called  sudden  ionospheric  disturb¬ 
ances  (SID),  such  as  short-wave  fade-out,  sudden  cosmic  noise  absorption, 
sudden  enhancement  of  atmospherics,  etc.,  are  manifestations  of  ioniza¬ 
tion  in  the  D- region  which  in  turn  is  associated  with  the  occurrence 
of  a  solar  flare.  Only  recently,  however,  has  the  nature  of  the  ionizing 
portion  of  the  solar  flare  radiation  come  to  light.  Measurements  (Friedman, 
Chubb,  Kupperian,  and  Lindsay  1958;  Chubb,  Friedman,  and  Kreplin,  1960a; 
Chubb,  and  others,  1960b;  Kreplin,  Chubb,  and  Friedman,  1962)  by  the 
Rocket  Astronomy  Group  at  the  l.S.  Naval  Research  Laboratory  ( NRL )  have 
indicated  that  this  radiation  in  the  hard  portion  of  the  spectrum  pene¬ 
trates  deeply  into  the  D-region  and  during  large  disturbances  apparently 
produces  ionization  at  altitudes  of  less  than  60  km. 

The  electron  concentration  is  influenced  not  only  by  the  intensity 
of  the  ionizing  radiation  but  also  by  the  recombination  rate  and  by 
the  electron  detachment  and  attachment  rates.  It  will  be  shown  later, 
however,  that  the  photodetachment  rate  is  sufficiently  large  so  that 
tin  (ffect  of  electron  attachment  to  molecular  oxygen  is  of  negligible 
importance  at  altitudes  above  75  km;  this  statement  is  true  during  the 
decay  phase  as  well  as  during  the  build-up. 

In  section  2  we  consider  the  NRL  measurements  of  solar  flare 
ladiation  and  their  implications  with  respect  to  the  D-region,  and  on 
the  basis  of  these  data  we  construct  a  model  of  the  electron  production 
prof  i  1° .  In  section  3  we  estimate  the  dissociative  recombination 
coefficient  and  construct  a  model  of  radio  signal  absorption  hased  on 
the  NRL  data. 
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The  Nature  and  Effects  of  Ionizing  Radiation  from  Solar  Flares 


2  . 


Since  World  War  II  several  theories  have  appeared  regarding  the 
nature  of  the  radiation  responsible  for  sudden  ionospheric  disturbances. 
Thr  f i i st  was  that  the  increase  in  ionization  is  caused  by  an  enhancement 
in  Lyman-a  emission  during  a  flare.  The  difficulties  inherent  in  this 
suggestion  are  now  apparent:  l,yman-a  would  have  to  increase  in  intensity 
by  a  factor  of  ~  10s  or  more  in  order  to  produce  the  high  frequency 
absorption  observed  during  sudden  ionospheric  disturbances.  The  detec¬ 
tion  of  X-rays  in  the  rocket  and  satellite  observations  of  solar  flares 
(Sunflare  1  and  II)  by  Friedman,  Chubb,  Kuppcrian,  and  Lindsay  (1938); 
Chubb,  Friedman,  and  Kreplin  (1960a);  Kreplin,  Chubb,  and  Friedman, 

(1962  ,  provides  a  basis  for  a  much  better  explanation  of  the  origin  oi 
sudden  ionospheric  disturbances.  The  X-ray  flux  is  spread  out  in  a 
spectrum  which  can  cause  sufficient  ionization  at  all  levels  of  the 
D-region  to  produce  the  observed  electron  concentrations.  It  has  also 
been  suggested  (Mitrn,  1960)  that  both  X-rays  and  Lyman-a  are  important 
in  producing  SID’s.  However,  lecent  satellite  measurements  of  the 
Lyman-a  Intensity  during  a  class  1+  solar  flare  ^Chubb  and  others,  1960b) 
indicate  that  no  such  enhancement  in  the  Lyman-a  intensity  occurs.  A 
larger  flare  may,  of  course,  produce  an  observable  Lyman-a  enhancement, 
but  it  is  doubtful  if  it  could  compete  with  X-radiation  in  its  ionospheric 
effects  at  any  level  of  the  D-region.  All  the  measurements  of  the  NRL 
group  Friedman,  Chubb,  Kuppcrian,  and  Lindsay,  1958;  Chubb,  Friedman, 
and  Kreplin,  1960a,  Chubb  and  others,  1960b)  indicate  that  the  electro¬ 
magnetic  solar  radiations  producing  the  enhanced  ionization  in  the  II- 
region  are  X-rays  of  wavelength  )  .  DA,  X-rays  of  wavelengths  ■  ->  8  A 
will  be  absorbed  almost  entirely  in  the  E-reglon  This  is  shown  in 
Table  111-1,  which  gives  the  intensity  at  several  wavelengths  required 
to  produce  given  electron  production  rates  at  two  levels  of  the  IJ-region 

A  Partial  X-ray  flux  spectrum  was  obtained  from  a  Sunflare  II 
rocket  Chubb,  Friedman,  and  Kreplin,  1960a)  launched  during  a  class  2+ 
flare  on  August  31,  1959  Infortunately ,  only  a  lower  limit  of  the 
intensity  in  the  2-8  A  range  was  observed  owing  to  saturation  of  some  of 
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Table  III-l 


INTENSITY  REQUIRED  TO  PRODUCE  GIVEN  ELECTRON  PRODUCTION 
RATES  AT  TWO  LEVELS  IN  THE  D-REGION 


Al ti tude , 

Electron 

Production 

Intensi ty 

,  erg  cm2  sec-1 

km 

Rate  q, 
cm-3  sec-1 

2  >  X  >  8  A 

X  '  8  A 

lq 

70  km 

30 

0.2 

^>1 

~10B 

90  km 

2500 

0.2 

-  1 

~10s 

the  sensors.  Because  of  this  gap  in  the  spectral  data  it  was  necessary 
to  estimate  the  missing  portion  by  extrapolating  from  other  data. 
Several  other  successful  rocket  flights  in  the  Sunflare  II  program 
yielded  essentially  complete  X-ray  spectra  in  the  range  2-60  A,  but 
because  they  were  made  during  a  small  flare  and  under  quiet  sun 
conditions,  the  cosmic  noise  absorption  data  were  inadequate  for  the 
construction  of  a  model  of  the  disturbed  D-region.  However,  since  we 
had  X-ray  spectra  for  several  lower  levels  of  X-ray  intensity,  it  was 
possible  to  extrapolate  the  shape  of  the  X-ray  spectrum  to  the  level 
of  X-ray  intensity  corresponding  to  the  class  2+  flare  considered  in 
this  report.  The  spectrum  so  constructed,  together  with  its  estimated 
upper  and  lower  limits  (the  broken  curves)  at  time  22:5-1  UT,  is  shown 
in  Figure  III-l,  in  addition,  counter  measurements  in  the  energy 
range  15  to  80  keV  were  also  made  at  two  different  times  during  the 
flight  These  data  indicated  a  softening  of  the  hard  X-radiation  as 
it  decreased  in  intensity  (Chubb,  Friedman,  and  Kreplin,  1960a). 


When  X-rays  are  absorbed  in  the  ionosphere,  they  produce  (a) 
primary  ionization  by  photoelectric  "stripping”  of  electrons  from 
neutral  molecules,  and  (b)  secondary  ionization  by  the  "stripped" 
electrons.  A  photon  of  wavelength  2  A  will  produce  about  200  free 
electrons  by  this  process.  The  electron  production  rate  can  be  written 
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HQ(X)  fi(X) 

dX  n  E 
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where  ^air(X)  is  the  absorption  cross  section  of  X-rays  in  air  as  a 
function  of  wavelength  X;  Q(X)  is  the  X-ray  photon  flux  density;  n  is 
the  neutral  particle  number  density;  E(X)  is  the  energy  of  a  photon  of 
wavelength  X;  and  Ee  is  the  mean  energy  required  to  remove  an  electron 
from  a  neutral  molecule  in  air  (~  32  eV).  Owing  to  absorption  by  the 
atmosphere,  the  X-ray  photon  flux  density  Q(X)  is  also  a  function  of 
altitude.  Figure  1 1 1 -2  shows  the  attenuation  of  I.-Q  and  X-radiation 


FIG.  111-2  ATTENUATION  OF  SOLAR  RADIATION  IN  THE  UPPER 
ATMOSPHERE.  THESE  CURVES  ARE  BASED  ON  THE 
ARDC  1959  MODEL  ATMOSPHERE  (MINZNER,  1959) 

of  several  wavelengths  in  the  upper  atmosphere.  Figure  III-3  shows  the 
electron- product ion-rate  profile  computed  from  the  data  presented  in 
Figure  1 1 1 -1 . 

We  shall  assume  in  our  model  that  the  electron  production  rate 
bears  the  same  functional  relationship  to  time  at  all  altitudes;  this 
is  equivalent  to  assuming  that  the  shape  of  the  X-ray  spectrum  is  time- 
independent  or  that  Q  can  be  separated  into  two  factors,  one  of  which 
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FIG.  MI-3  ELECTRON  PRODUCTION  RATES  AT  VARIOUS  ALTITUDES 
AS  COMPUTED  FROM  THE  SOLID  CURVE  IN  FIG.  Ill-l 

depends  only  on  wavelength  X  and  altitude  h,  and  the  other  only  on 
time 

Q(X,  h,  t)  =  A(X,  h)T(t)  (2) 

We  can  establish  the  start  of  the  X-ray  emission  at  the  time  of  the 
commencement  of  sudden  cosmic  noise  absorption  (see  Figure  III-4)  and 
the  maximum  of  the  X-ray  emission  curve  at  roughly  the  time  of  minimum 
cosmic  noise  intensity.  This  is  an  oversimplification,  since  there  is 
no  reason  to  believe  that  the  entire  X-ray  spectrum  bears  exactly  the 
same  functional  relationship  to  time. 
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Gala  reproduced  through  courtesy  of  U  S  Naval  Research  Laboratory 

FIG.  1 1 1-4  INTENSITY  AT  LOS  ANGELES  OF  18-Mc  COSMIC  NOISE  DURING 
THE  OCCURRENCE  OF  A  CLASS  2*  SOLAR  FLARE  ON 
31  AUGUST  1959.  THE  RECORD  INDICATES  THE  EXTRAPOLATION 
TO  ITS  ORIGIN  OF  THE  LARGER  OF  THE  SUDDEN  COSMIC 
NOISE-ABSORPTIONS 


3 .  Model  of  Ionization,  Recombination ,  and 
Radio  Signal  Absorption 


A  sudden  ionospheric  disturbance  may  be  divided  into  two  phases; 
the  build-up  phase  and  the  decay  phase.  The  former  is  characterized  by 
a  rapid  increase  in  electron  density  due  to  the  increasing  solar  X-ray 
photon  flux  and  thus  the  electron  production  rate;  the  latter,  on  the 
other  hand,  is  characterized  by  a  gradual  deciease  in  electron  density. 
The  electron  density  maxima  do  not  necessarily  occur  simultaneously  at 
all  altitudes,  however,  since  there  is  at  present  no  known  reason  for 
the  intensity  maxima  of  ail  parts  of  the  the  X-ray  spectrum  to  occur 
simultaneously.  In  spite  of  this,  we  shall  assume  in  our  model  that 
the  electron  densities  do  simultaneously  reach  their  maxima  at  all 
altitudes  in  the  I)- region  and  that  each  maximum  occurs  at  the  same  time 
as  the  minimum  in  the  cosmic  noise  intensity  curve  (Figure  1 1 1  - 4 )  .  The 
time  rate  of  change  of  electron  density  is  related  to  the  production 
and  removal  processes  by  the  well-known  equation 
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in  which  Ng  is  the  electron  density;  q  is  the  electron  pioduction  rate; 
>  is  the  negative  ion  to  electron-density  ratio  (X  =  N~/N  );  is  the 
dissociative  recombination  coefficient;  and  Q.  is  the  mutual  ionic 
neutralization  coefficient.  The  meaning  of  the  symbol  X  (wavelength 
or  negative- ion  to  electron  concentration  ratio)  in  the  following  will 
be  clear  from  the  context. 


The  differential  equation  obeyed  by  the  negative  ion  to  electron 
density  ratio  X  is  (Bates  and  Massey,  1946) 


1  + 
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where  ft  is  the  electron  attachment  coefficient;  ,  is  the  photodetachment 
rate;  y  is  the  electron  collisional  detachment  coefficient;  n  is  the 
neutral  particle  density;  and  m  is  the  molecular  oxygen  concentration 
(taken  to  be  1/5  n  in  our  computations).  In  the  model  to  be  developed 
in  this  report  the  neutral  particle  concentrations  given  in  the  A  Ft  DC  1959 
model  atmosphere  (Minzner,  Champion,  and  Pond,  1959)  are  adopted. 

The  second  term  in  brackets  on  the  right-hand  side  of  equation  (1 
is  of  importance  only  in  the  lower  D-region.  Recent  estimates  of  the 
collisional  detachment  coefficient  y  based  on  polar  cap  absorption 
(Bailey  and  Branscomb,  I960)  an  of  the  order  y  2  x  i()“17  cm3,  sec 
However,  other  workers  (Phelps  and  Pack,  1961)  ban  Found  evidence 
that  y  s-  10” 1 7  cm 3/sec  for  02~  ,  and  that  N2  is  an  ineffect. ve  detaching 
agent.  Whitten  and  Poppoff  (1962  have  suggested  that  an  associative 
detachment  mechanism  may  b<  dominant  Bo  that  as  it  may,  our  investiga¬ 
tion  showed  that  collisional  detachment  is  not  very  important  in  the 
solar  flare-disturbed  D- region.  Tin  third  and  Fourth  terms  in  brackets 
are,  according  to  our  computations,  approximately  equal  in  magnitude 
but  of  opposite  sign;  hence  they  tend  to  cancel  The  negative  sign  of 
the  third  term  in  brackets  N<_  -  vyj  is  due  to  the  Fact  that 

I)  \  Crain>  1961  Because  ol  its  verv  slight  dependence  on  N 
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and  q,  X  is  a  slowly  varying  function  of  time1  and  can  be  approximated 
by 

.  pm 

x  35  —  (4a) 

Because  of  the  weak  time-dependence  of  X,  the  last  term  on  the  right-hand 
side  of  equation  3  can  be  neglected.  The  most  recent  investigations 
of  electron-02  attachment  (Chanin,  Phelps,  and  Biondi  1959;  Holt,  1959) 
indicate  that  in  the  D-regicn  it  is  primarily  a  three-body  process  and 
that  the  coefficient  K  where 


P  =  Km 


(5) 


has  the  value  K  =  (2  ±  lj  x  10-30 
250°K .  Equation  (4a)  can  then  be 


cm6/ sec 
written 


lor  an  electron  temperature  of 


The  photodetachment  rate  p  was  computed  using  recent  measurements 
(Burch,  Smith,  and  Branscomb  1958)  of  photodetachment  cross  section  for 
02  and  solar  radiation  flux  measurements  (Kuiper,  1953).  These  yield 
a  value  of  p  =  0.35  see-* .  The  value  of  the  electron  affinity  of  02 
derived  by  Burch  and  others  from  the  cross  section  data  was  0.15  eV 
compared  with  the  value  of  0.46  eV  obtained  by  Phelps  and  Pack  (1961). 
Since  the  photodetachment  cross  sections  are  expected  to  be  decreasing 
functions  of  the  threshold  energy  of  the  process,  our  computed  value  of 
t  may  be  too  large.  Hence  X  was  also  computed  for  values  of  p  of 
0.10  sec  1  and  0.035  sec-1;  the  x  profiles  are  shown  in  Figure  1II-5. 

The  value  >  =  0.1  sec-1  is  provisionally  suggested  as  the  most  probable 
value  (see  section  IVA  i  . 


The  insensitivity  of  X  to  changes  in  q  and  N  indicates  that 
to  equilibrium  of  the  attachment  and  detachment  processes  is 
This  is  merely  the  result  of  relatively  large  attachment  and 
coefficients. 


the  approach 
quite  rapid, 
detachment 
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FIG.  Ill- 5  NEGATIVE  ION-ELECTRON  CONCENTRATION  RATIO  PROFILES 
FOR  THREE  POSSIBLE  VALUES  OF  THE  PHOTODE  UCHMENT 
RATE 


If  the  values  of  K,  K ,  j,  i^,  a^ ,  the  neutral  particle  concent t  i- 
tion  profile,  and  the  functional  relationship  to  altitude  and  time  of 
the  electron  source  function  q(h,t)  are  known,  it  is  possible  to  eomput, 
the  electron  concentration  at  any  altitude  h  by  means  of  the  approximat. 
form  of  equation  (3): 


Using  the  foregoing  as  a  basis,  we  shall  now  construct  a  model  of 
the  D-region  during  sudden  ionospheric  disturbances.  To  facilitate  this 
and  to  obtain  an  estimate  of  the  dissociative  recombination  coefficient 
we  shall  utilize  the  N’RL  Sunflare  II  data  and  sudden  cosmic  noise  absorp¬ 
tion  data  presented  m  Figure  III--!.  In  essence  our  task  is  to  find  a 
functional  form  for  q  and  a  value  of  the  dissociative  reeombi nat ion 
coefficient  which,  when  inserted  in  equation  (3a),  vield  electron 
concentration  profiles  at  various  times  during  the  build-up  phase  that 
agree  with  the  observed  reduction  in  18  Me  cosmic  noise  intensity.  The 
mutual  ionic  neutralization  coefficient  Cl  does  not  enter  the  computation 
because  it  is  probably  smaller  than  CL  (Crain,  1961)  and  because  of 
the  fact  that  ttie  negative  ion-electron  concentration  ratio  at  most  is 
of  order  unity  at  altitudes  of  interest 


The  relative  cosmic  noise  intensity  is  given  by  the  equation 


Id 

—  =  cxPi-2  (>  (kd  -  kn)  dhj  (6) 

n 

where  1^  and  I  are  the  18  Me  cosmic  noise  intensities  corresponding  to 
the  disturbed  and  normal  D-region,  respectively,  and  /  indicates 
integration  over  that  segment  of  the  ray  path  which  lies  in  the  D-region. 
The  absorption  coefficient  k  foi  the  case  in  which  a  »  v  is  given 
approximately  by 


k 


(7) 


Mol  mud ,  1959),  neglecting  gyromagnetie  splitting  which  is  of  little 
importance  at  18  Me;  0  is  the  plasma  frequency,  j.  ,  the  signal  frequency, 
\ ,  the  elec t ron— neu t ral  particle  collision  frequency,  and  c  the  velocity 
of  light  in  free  space. 


In  constructing  this  model  we  considered  electron  collision 
frequencies  suggested  by  Phelps  and  Pack  (1959),  by  Nicolet  (1959),  by 
Kane  (i960),  and  by  Barrington  and  Thrane  (1962).  Those  proposed  by 
Phelps  and  Pack  were  based  on  drift  velocity  measurements  on  electrons 


in  Na  and  probably  constitute  a  lower  limit,  while  those  due  to  Nicolet 
were  based  on  the  results  from  a  microwave  interaction  experiment  carried 
out  by  Anderson  and  Goldstein  (1956).  The  latter  measurements  undoubtedly 
corresponded  to  much  larger  electron  temperatures  than  was  realized  at 
the  time  (Forma to  and  Gilardini,  1959)  and  thus  yielded  collision 
frequencies  which  are  too  large.  Lying  between  these  models  are  the 
experimental  results  of  Kane  (1960)  which  were  based  on  in  situ  measure¬ 
ments  of  high  frequency  radio  wave  absorption,  and  of  the  cross  modula¬ 
tion  results  of  Barrington  and  Thrane  (1962).  Since  these  results  agree 
quite  well  for  a  quiet  ionosphere,  we  employed  them  in  the  computation 
of  the  dissociative  recombination  coefficient. 

The  X-ray  spectrum  given  in  Figure  III-l,  together  with  equations (l) 
and  (2),  may  be  used  to  compute  the  electron  production  rate  profile 
as  a  function  of  time  if  we  assume  a  model  for  X-ray  intensity  time- 
dependence  T(t) . 

Several  simple  functional  forms  of  T(t),  as  well  as  several  values 
of  the  dissociative  recombination  coefficient  C^,  were  used  in 
equation  (3a),  which  was  then  integrated  numerically  in  order  to  obtain 
the  electron  concentration  profiles  at  various  times  during  the  build¬ 
up  phase.  These  profiles  were  then  used  in  conjunction  with  equations  (6) 

and  (7  to  compute  the  relative  cosmic  noise  intensity  I  I  as  a 

d  n 

function  of  time.  Upon  comparson  of  I ^ ( t )/ I  with  the  sudden  cosmic 

d  n 

noise  absorption  data  presented  in  Figure  1 1 1-1 ,  it  was  immediately 
evident  that  only  one  of  the  trial  analytical  forms,  i.e.,  T(t)  =  at™, 
was  a  reasonably  valid  choice.  After  various  values  of  m  in  the  range 
i  £  m  <  1  were  tried,  it  was  found  that  m  =  1/2  best  agreed  with  the 
sudden  cosmic  noise  absorption  data.  The  value  of  j  that  best  agreed 
with  the  SCNA  data  was  =  (2^  x  10-7  cm3/sec ;  the  uncertainty  in  a 
reflects  the  uncertainties  in  the  X-ray  spectrum  and  the  electron 
collision  frequency  profile. 
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The  decay  phase  was  then  treated  in  a  similar  manner.  It  was  found 
that  the  functional  form  T(t)  =  a'(t0  -  t):  where  t0  =  11.7  min,  and 
a  value  of  =  2  x  10-7  cm3/ sec  were,  when  used  to  compute  the  relative 
cosmic  noise  intensity  I^tj/I^  in  the  manner  outlined  previously,  in 
reasonable  agreement  with  the  sudden  cosmic  noise  absorption  observa¬ 
tions  until  15  minutes  after  commencement  of  the  flare.  The  solid  curve 
in  Figure  III —6  shows  the  analytical  form  for  T  t)  used  in  these  conputatiais. 


FIG.  111-6  MODEL  OF  SOLAR  RADIATION  INTENSITY  0  -  2  to  8A)  AS  A 
FUNCTION  OF  TIME  FOR  THE  FLARE  OCCURRING  ON 
31  AUGUST  1959.  THE  SOLID  LINE  REPRESENTS  THE 
FUNCTION  ot'j,  0  _  »  _  5  MIN.  AND  o'  (11.7  t)2,  t  s  5  MIN. 

THE  BROKEN  LINE  REPRESENTS  THE  ESTIMATED  X-RAY 
INTENSITY 


The  18  Me  absorption  exponent  2  f  (k,  -  k  dh  is  shown  as  a 

'  L)  '  d  n 

function  of  time  in  Figure  III —7 .  The  solid  curve  was  obtained  by 
substituting  into  equation  (3a)  the  analytical  lorm  of  T(t)  and  the 
values  of  previously  derived,  then  numerically  integrating  the 
equation  to  obtain  the  electron  concentration  profile,  and  substituting 
the  electron  concentrations  N^(h,  t)  so  obtained  into  equation  (7). 
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FIG.  111-7  18-Mc  COSMIC-NOISE-ABSORPTION  EXPONENT  2  kdh  DURING  THE 
BUILD-UP  AND  DECAY  PHASES  OF  THE  FLARE  OCCURRING  AT 
22.50  UT  ON  31  AUGUST  1959.  THE  BROKEN  CURVE  REPRESENTS 
THE  OBSERVED  ABSORPTION  AND  THE  SOLID  CURVE 
REPRESENTS  THE  COMPUTED  ABSORPTION  USING  THE  FUNCTIONS 
at'2  (0  <  t  <  5  min)  AND  a  '  (1 1.7  “  t)2  (5  min  <  t  <  117  min)  TO 
REPRESENT  THE  ELECTRON  PRODUCTION  RATE.  THE  SYMBOL  k 
IS  EQUAL  TO  k  .  -k 

o  n 

The  profiles  of  [k^(t)  -  j  were  then  integrated  over  the  ray  path 
through  the  D-region  to  obtain  the  absorption  exponent.  It  is  evident 
from  Figure  III —7  that  the  functional  relationship  to  time  of  X-ray 
intensity  T(t)  actually  increased  faster  than  at*/2  during  the  initial 
portion  of  the  build-up  phase  and  decreased  faster  than  a'(t0  -  t)z  during 
the  decay  phase.  This  is  indicated  qualitatively  by  the  broken  curve 
in  Figure  III-6;  the  broken  curve  agrees  qualitatively  with  X-ray 
intensity  time-variation  data  obtained  for  two  different  flares  by 
satellite  1960  eta  2{R.  W .  Kreplin,  private  communication).  We  wish 
to  reiterate  that  the  factorization  of  Q(\,  h,  t)  in  accordance  with 
equation  (2)  is  quite  artificial.  It  is  known,  for  example,  that  the 
X-ray  sprectral  form  changes  with  time  during  the  decay  phase  of  a 
flare;  this  undoubtedly  accounts  for  part  of  the  disagreement  between 
the  observed  and  computed  absorption  exponents. 


22 


t*'r»AT'C« 


Figure  III —8  shows  the  electron  concentration  profile  at  22:54  UT 
as  computed  by  the  method  outlined  previously.  The  corresponding  absorp¬ 
tion  coefficient  obtained  by  substituting  the  data  contained  in 
Figure  III-8  into  equation  (7)  is  presented  in  Figure  III-9.  This 
absorption  coefficient  profile  is  in  rough  agreement  with  the  model 
proposed  by  Mitra  (1960);  the  maximum  at  about  80  km  occurs  about  10  km 
higher  than  is  proposed  in  the  high  frequency  absorption  model  of 
Kanellnkos  and  Villard  (1962). 

B,  Ionization  by  Background  Solar  X-rays 

The  significance  to  normal  D-region  ioni/atlon  of  recently  reported 
values  of  background  solar  X-ray  fluxes  (Kreplln,  1961)  has  largely  been 
ignored.  The  purpose  of  this  section  is  to  present  the  results  of  a 
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FIG.  1 1 1 -9  18  Me  ABSORPTION  COEFFICIENT 
PROFILE  CORRESPONDING 
TO  THE  ELECTRON  CONCENTRATION 
PROFILE  PRESENTED  IN  FIG.  1 11-7 


FIG.  111-8  COMPUTED  ELECTRON 

CONCENTRATION  PROFILE 
AT  TIME  22:54  UT 


1961  study  to  assess  the  possible  role  of  background  solar  radiation 
in  the  formation  of  the  D-region. 

Current  concepts  of  normal  D-region  ionization  processes  seem  to 
favor  the  model  proposed  by  Nicolet  and  Aikin  (1960),  although  other 
workers  (Inn,  1961a,  b;  Bourdeau,  1961)  have  offered  feasible  alternative 
hypotheses.  There  seems  to  be  little  argument  about  the  dominant  role 
of  eosmie  radiation  in  the  production  of  lower  D-region  ionization.  The 
crux  of  the  present  controversy  is  the  lack  of  verified  data  on  the 
relative  influence  of  solar  Lyman-Ct  radiation,  nitric  oxide,  molecular 
oxygen,  cosmic  radiation,  particulate  matter,  and  X-radiation  in  producing 
upper  D-region  ionization.  Nicolet  and  Aikin  (1960)  concluded  that  free- 
eleetron  concentrations  between  approximately  80  and  85  km  are  produced 
principally  as  a  result  of  the  ionization  of  nitric  oxide  by  solar 
Lyman-Q!  radiation.  This  conclusion  is  the  result  of  calculations  employ¬ 
ing  a  theoretical  vertical  distribution  of  nitric  oxide,  an  assumed 
dissociative  nitric  oxide  ion-electron  recombination  coefficient,  and 
the  solar  radiation  measurements  available  before  1960.  Recent  work 
by  Gunton  and  Inn  (1961)  and  by  the  authors  (section  IV  C,  this  report) 
shows  that  the  dissociative  recombination  coefficient  assumed  by  Nicolet 
and  Aikin  may  have  been  too  small  by  two  to  three  orders  of  magnitude. 

This  means  either  that  the  nitric  oxide  concentrations  are  higher  than 
those  calculated  by  Nicolet  and  Aikin  or  that  another  ionization  mecha¬ 
nism  must  be  found. 

Inn  (1961a, b)  suggests  a  mechanism  based  on  the  ionization  of  vibra¬ 
tional  ly  excited  ground-state  oxygen  molecules  by  solar  Lyman-a.  Bourdeau 
(1961)  suggests  that  cosmic  ray  ionization  modified  by  diffusion  of  ions 
to  particulate  material  in  the  60  to  80  km  region  can  account  for  the 
observed  conductivity  profiles. 

We  should  like  to  suggest  another  possibility  by  considering  the 
X-ray  spectra  reported  by  Kreplin  (1961)  as  being  characteristic  of 
background  solar  radiation,  at  least  during  periods  of  sunspot  maxima. 
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Kreplin  reports  X-ray  fluxes  of  the  order  of  10"3  erg/cm2  sec  in  the  2-8  A 
range  during  quiescent  solar  conditions.  These  values  are  considerably 
higher  (by  two  or  three  orders  of  magnitude)  than  those  utilized  for 
the  quiet  sun  condition  by  Nicolet  and  Aikin;  therefore,  an  evaluation 
of  the  possible  effect  of  this  level  of  solar  radiation  on  the  I, -region 
seems  appropriate.  This  suggestion  was  made  years  ago  by  Muller  (1935), 
again  by  Rawer  (1952),  and  more  recently  by  Chamberlain  (1961).  Until 
now,  however,  data  have  not  been  available  to  substantiate  the  suggestion. 

The  data  reported  by  Kreplin  (1961)  from  flight  NN  8.69  CK  included 
values  of  solai  X-radiation  for  three  energy  ranges  and  were,  therefore, 
ustd  to  prepare  an  estimated  quiet  sun  X-ray  spectrum  as  shown  in 
Figure  I I I — 1 0 .  The  data  from  flight  NN  8.68  CF ,  though  not  as  complete, 
substantiate  the  data  from  flight  NN  8.69  CF . 


FIG.  Ill- 10  BACKGROUND  SOLAR  X-RAY  SPECTRUM 
CONSTRUCTED  FROM  DATA  OBTAINED 
BY  SUNFLARE  II  PROBE  NN  8.69  CF 
AT  15:59  UT  ON  8  14  59 
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Measurements  of  background  X-ray  spectra  and  intensities  at  sunspot 
maximum  reported  by  other  observers  tend  to  confirm  the  NRL  data  and 
our  estimated  spectrum.  The  spectrum  measured  by  the  SL  40  rocket 
(9.27.61)  is  in  close  agreement  with  our  estimate  in  the  4-9  A  region 
(Pounds,  1962).  The  measurement  of  intensities  by  Pounds  and  Sanford 
(1962)  during  the  period  1959  to  1961  varies  from  10"3  ergs  cm-2  sec-1 
in  the  4-10  A  range  to  1.5  x  10" 3  ergs  cm"2  sec"1  in  the  2-8  A  range. 
Mandel'shtam,  et_  a_l .  ,  (1962)  report  rocket  measurements  of  1.6  x  10“4 
ergs  cm  sec  1  in  the  0-8  A  range  to  7.3  x  10“4  ergs  cm"2  sec-1  in  the 
0-10  A  range  in  1959.  Reported  measurements  of  0-10  A  background  X- 
radiation  by  satellite  observatories  in  succeeding  years  vary  from  less 
than  0.6  x  10"3  ergs  cm"2  sec"1  (Kreplin,  Chubb,  and  Friedman,  1962) 
to  1.5  x  10-4  ergs  cm"2  sec"1  (Neupert,  et^  al. ,  1962).  Note  that  all 
of  these  measurements  are  within  the  range  found  to  be  significant  for 
D-region  ionization  by  the  calculations  described  in  this  section. 

Using  the  spectrum  of  Figure  III-10,  the  appropriate  solar  zenith 
angle  (45  ),  the  appropriate  X-ray  cross  sections,  and  ARDC  1959  model 
atmosphere,  a  profile  of  electron  production  rates  was  calculated.  It 
is  shown  in  Figure  III-ll. 


Since  the  radiation  measurements  were  made  during  quiet  solar  con¬ 
ditions,  it  is  a  good  approximation  to  consider  the  D-region  to  have 
been  in  a  steady  state.  The  electron  density  N  (h)  at  altitude  h  is 
then  given  by 


Vh)  = 


-iQi) 


l  * 


(i  +  X)(a0  +  xat ) J 


(i) 


in  which  is  the  negative  ion  to  electron  concentration  ratio,  Ct^  is 

the  dissociative  recombination  coefficient,  and  a  is  the  mutual  ionic 

neutralization  coefficient  It  can  be  shown  that,  i f  we  restrict  the 

discussion  to  altitudes  above  70  km,  omission  of  XCt  will  introduce 

an  error  of  less  than  half  an  order  of  magnitude  in  N  (h)  at  the  lower 

e 

altitudes,  thus  simplifying  equation  (1).  From  our  assumption  that 
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FIG.  111-11  ELECTRON  PRODUCTION  RATE  AND  CONCENTRATION 
PROFILES  CORRESPONDING  TO  THE  X-RAY 
SPECTRUM  SHOWN  IN  FIG.  111-10 

the  ionization  is  caused  by  solar  X-rays,  it  follows  that  the  ionic 
processes  would  be  similar  to  those  in  the  same  region  produced  by  solar 
flare  radiation.  We  considered  this  problem  in  an  earlier  paper  (Whitten 
and  Poppoff,  1961),  deducing  a  value  of  =  3  x  10"7  cm6/ sec  and  a  X 
profile  given  by  the  curve  labeled  o  =  0.10  in  Fipure  III-5. 

Using  q (h )  from  Fipure  III-ll,  the  X  profile  referred  to  above, 
and  the  above  value  of  a  ,  the  electron  concentration  profile  of 
Fipure  III-ll  is  obtained.  The  values  obtained  are  within  the  spread 
of  observed  values  for  D-region  electron  concentrations  at  these  altitudes 
Ratcliffe,  I960).  The  X-ray  intensity  could  be  an  order  of  magnitude 
lower  and  still  he  adequate  to  account  for  some  of  the  lower  estimates 
of  D-region  electron  concentrations.  If  the  spectrum  shown  in-Figure  1 1 1 -10 
is  'softened'  so  that  the  X-ray  intensity  at  X  =  4  A  is  about  10 ~s  erp, 
cm2  sec  A,  the  100  cm  3  level  in  the  electron-density  profile  (Figure  III-ll) 
is  raised  to  about  76  km. 
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The  electron  concentration  profile  of  Figure  III-ll  is  not  intended 
to  serve  as  a  model  for  the  D-region.  It  does,  however,  illustrate  the 
possible  role  of  X-radiation  in  the  altitude  range  70  to  90  km.  It 
also  shows  that,  if  the  solar  radiation  flux  measurements  used  for  the 
calculations  arc  representative  (to  within  an  order  of  magnitude), 
electron  concentrations  in  the  upper  D-region  can  be  accounted  for  with¬ 
out  recourse  to  hypotheses  requiring  Lyman-a  ionization.  Preliminary 
calculations  indicate  that,  if  the  electron  concentrations  estimated 
above  are  combined  with  estimates  of  ionization  produced  by  cosmic  radia¬ 
tion,  the  classical  D-region  electron  profile  can  be  obtained  which 
includes  a  peak  or  inflection  (depending  on  relative  values  of  cosmic 
and  X-ray  fluxes)  at  middle  D-region  altitudes.  A  reasonable  model 
must,  of  course,  include  all  valid  hypotheses  of  D-region  electron 
production  processes.  The  inclusion  of  Lyman-a  ionization  processes, 
as  well  as  cosmic  ray  and  X-ray  ionization  processes,  would  probably 
result  in  electron  profiles  with  several  inflections  or  peaks  (or  at 
least  changes  of  gradient)  whose  position  would  vary  with  diurnal, 
seasonal,  and  solar-cyclical  patterns. 

C .  Ionization  by  Auroral  Electrons 

1.  Introduction 


For  many  years,  observations  have  been  made  of  "soft 


Lift.  U  W" 


magnetic  radiation  at  low  altitudes  in  the  auroral  zone.  Van  Allen 
,1955)  suggested  that  the  observed  radiation  was  bremsstrahlung  generated 
by  low  energy  electrons.  Following  this  suggestion,  Chapman  and  little 
1907  and  Chapman  1959)  prepared  models  of  the  lower  ionosphere  which 
showed  that  this  interpretation  was  consistent  with  observations  of 
mdio  absorption  in  the  auroral  zone.  Meanwhile,  observations  were  made 
by  Anderson  and  Kncmark  (i960)  of  the  bremsstrahlung  spectrum  which 
furnished  the  basis  for  estimates  of  the  spectrum  of  the  generating 
electrons.  Although  the  source  of  the  electrons  is  still  somewhat  con¬ 
troversial,  it  has  been  shown  that  qualitative  correlations  exist  with 
geomagnetic  activity  and  with  solar  activity. 
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Close  correlations  with  cosmic  noise  absorption  have  been  observed. 
This,  of  course,  suggests  a  possible  approach  to  the  study  of  atmospheric 
relaxation  processes,  similar  to  the  SID  study  made  earlier  in  this 
project  and  described  in  III  A  (Whitten  and  Poppoff,  1961).  However, 
the  area  of  electron  bombardment  is  relatively  small  and  not  stationary. 
Thus  simultaneous  observations  of  the  radiation  from  balloons  and  of 
cosmic  noise  absorption  by  ground-based  riometers  were  very  difficult 
to  obtain  and,  in  fact,  did  not  exist.  Recently,  however,  Brown  and 
Campbell  v1962)  reported  simultaneous  observations  of  radio  noise 
absorption  and  radiation  flux  during  a  magnetic  disturbance  at  College, 
Alaska  on  June  25,  1961.  The  existence  of  these  observed  data  provided 
the  opportunity  to  attempt  the  construction  of  a  model  of  the  ioniza¬ 
tion  produced  by  auroral  electrons  and  the  resulting  electron  and  ion 
profiles.  Although  the  model  could  not  be  completed  for  this  report, 
a  considerable  portion  of  the  work  has  been  done  and  will  be  described 
in  this  section. 


2.  Ionization  by  Electrons 


The  procedure  followed  was  generally  that  suggested  by  Chamberlain 
1,1961 /(  which  is  essentially  to  consider  the  incident  electrons  as  a 
plane  source  normally  incident  to  the  atmosphere.  By  making  such  an 
assumption,  the  extensive  work  of  Spencer  (1959)  at  the  National  Bureau 
of  Standards  could  be  utilized.  The  weakness  of  this  assumption  would 
seem  to  lie  in  the  obvious  fact  that  the  electrons  are  not  normally 
incident  but  arrive  with  an  angle  of  incidence  controlled  by  the  geo¬ 
magnetic  field  in  the  bombardment  area.  Quite  obviously,  a  calculation 
that  considered  these  factors  would  be  beyond  the  resources  of  this 
study.  However,  calculations  reported  by  Pedersen  (1962)  show  that  the 
penetration  altitude  would  not  change  appreciably  with  angie  or  inci¬ 
dence.  For  example,  he  notes  that  100  keV  electrons  at  vertical  inci¬ 
dence  would  penetrate  to  80  km,  whereas  at  angles  of  incidence  30°  and 
60  from  the  vertical,  they  would  penetrate  to  81  and  85  km,  respectively 

Brown  and  Campbell  (1962)  reported  an  electron  flux  of  10s  per 
cm3  sec  in  the  energy  range  greater  than  50  keV  at  the  peak  of  the 
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ionization  event.  In  making  their  estimate  an  electron  spectrum  of 
Ne  =  k  E-5  was  assumed  and  the  computing  procedure  reported  by  Anderson 
and  Enemark  (i960)  was  utilized  to  obtain  electron  flux  estimates  from 
measured  X-ray  fluxes.  An  area  of  bombardment  of  15,000  to  -40,000  km2 
was  estimated. 

I'si  ng  the  above  data,  an  extension  of  the  National  Bureau  of 
Standards  work,  and  the  AROC  1959  model  atmosphere,  an  altitude  profile 
of  the  peak  ionization  rate  was  computed.  The  energy  range  was  extended 
to  10  keV  in  the  low  energy  and  to  200  keV  in  the  high  energy  end  of  the 
reported  spectrum.  The  results  are  illustrated  in  Figure  III-12  and  III-13. 


FIG.  111-12  IONIZATION  RATE  vs.  ALTITUDE  (PER  ENERGY  INTERVAL) 
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FIG.  111-13  TOTAL  IONIZATION  RATE  v5.  ALTITUDE  FOR  LOW  ENERGY 
CUTOFF  AT  10  kev  AND  20  kcv 


Study  of  the  calculated  curves  immediately  reveals  an  obvious 
difficulty.  Because  oi  the  steep  energy  spectrum  of  the  incident 
electrons,  the  lowest  energy  portion  is  always  dominant  in  determining 
the  peak  in  the  ionization  rate  profile.  Therefore,  the  low  energy 
cutoff  point  should  be  considered  very  carefully  in  the  preparation  of 
an  adequate  model . 

Although  electrons  have  been  detected  with  energies  as  low  as  8  keV 
(Davis,  Berg,  and  Meredith,  lf)60j,  recent  analysis  of  auroral  bremsst rahlung 
measurements  by  Winckler,  Bhavsar,  and  Anderson  1962)  suggests  that  the 
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electron  energy  is  concentrated  in  the  20-100  keV  region,  with  a 
maximum  intensity  in  the  50-80  keV  region  of  the  electron  energy 
spectrum. 

The  crux  of  the  problem  for  this  study  would,  of  course,  be  the 
absorption  of  radio  frequency  energy  by  the  D-region  The  maximum  in  the 
ionization  rate  profile  not  only  increases  in  magnitude  with  decreasing 
elic.tron  tnttgy,  but  also  in  altitude.  Since  collision  frequencies 
decrease  with  altitude,  a  compensating  decrease  should  also  occur  in 
radio  absorption  coefficients.  Thus,  the  effect  of  imprecise  knowledge 
of  the  low  energy  cutoff  may  not  bo  critical  Reid  (1961),  in  studying 
proton  ionization  effects  during  polar  radio  blackout  events,  found  that 
lack  of  knowledge  of  low  energy  proton  cutoff  points  was  not  very 
i mportant 

3  Ionization  by  Brcmsst rahl ung 

Obviously  (see  Figure  III-12)  the  altitudes  at  which  incident 
primary  electrons  are  absorbed  and  bremsstrahlung  is  produced  va  -y 
greatiy  with  energy,  from  about  75  to  100  km  for  the  energy  range 
considered.  Inasmuch  as  the  penetration  altitude  of  the  bremsstrahlung 
is  a  function  of  the  source  altitude,  this  should  be  considered  in  the 
calculation  of  bremsstrahlung  ionization  rate  profiles  In  these 
calculations,  it  was  assumed  that  bremsstrahlung  is  produced  at  the 
al t i tude  of  peak  absorption  of  the  incident  primary  electrons  in  the 
energy  ranges  used  in  preparing  the  curves  of  Figure  1 1 1-12. 

Profiles  of  ionization  produced  by  bremsstrahlung  g  nerated 
by  electrons  in  each  of  the  several  election  energy  ranges  were  computed. 
The  calculated  ionization  rates  were  then  summed  for  each  altitude  and 
a  Profile  was  prepared  of  ionization  rates  produced  by  bremsstrahlung 
generated  by  the  entire  spectrum  of  electrons. 

Bremsstrahlung  spectra  were  computed  by  following  the  treatment 
described  by  Anderson  and  Enema.*  (I960)  which  can  be  summarized  us 
foil  own  . 
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The  number  of  photons  produced  by  an  electron  of  kinetic  energy 

E  in  passing  through  a  thickness  dg  g  cm"2  of  material  having  radiation 
length  L  g  cm-2  is 

K 


d' 


dvn  =  f- 


v="  /h 


dv 


R  .  =e  ,  h 

a 


a 


V  =  frequency  of  photon 

Ea  =  chosen  discrimination  level  or  the  lower 
limit  of  energy  interval  considered. 


The  number  of  photons  with  energies  above  the  discrimination  level  K 
produced  by  an  electron  initially  of  energy  E  traversing  its  entire  ‘ 
range  is 


n(Ea;  E)  =  /  d  a  (?)  =  f-  /  In  L.  d. 

c  LR  Ea 

J  a 

Ising  the  approximate  range-energy  relation  R  =  K  2000,  where  K 

is  expressed  in  keV  and  R  in  g  cm"2  ,  the  above  equation  can  be  integrated 
to  give 

n  h  2000  L  ‘ E  ln  E  E  -  1  '  E  J 

It  a 

The  expression  for  another,  higher  discrimination  level,  E  would 
be  identical  except  for  the  subscripts  The  number  of  photoL  in  the 
energy  interval  Ea  to  Eq+j  produced  by  N  electrons  of  energy  E  is 


n(E 


E 

a+1 


*tE) 
2000  I 


K  ln 


"  (E.. 


ai  1 


-  Ka^j 


The  number  of  photons  produced  in  each  of  a  series  of  energy  intervals 
Iron,  1-2  kcV  up  to  the  energy  interval  of  the  primary  electrons  was 
the,,  computed  for  each  electron  energy  interval  used  in  Figure  III-12 
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The  bremsstrahlung  flux  at  the  distance  R  from  a  large  isotropic 
radiation  can  be  computed  as  follows;  (see  Figure  III-14).  The  eontri 
bution  from  each  unit  increment  is 


R2 

where  i0  is  the  number  of  photons  emitted  per  unit  source  area  and  r 
is  the  distance  from  the  center  of  the  source  to  the  unit  source.  A 
and  B  are  attenuation  and  scattering  terms  respectively. 


UNIT  INCREMENT 


/ 


/ 

-  UNIT 
VOLUME 
RECEPTOR 

OA  3624-33 

FIG.  III. 14  SCHEMATIC  OF  RADIATION  GEOMETRY 
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The  attenuation  term, 


A  =  exp  (  -  /  pT(E)dR), 
o 


can,  for  this  application,  be  expressed  as 


A  =  exp  (-  —  M  J) 
Po  d 


The  scattering  term  is  approximated  as 

M_ (E)  R 


B  =  1  + 


Pa(E)  ;  mt 


/  M_(E)dR 


and  can,  for  this  case,  be  expressed  as 


M  ti 

B  =  1  +  —  —  M  £ 

Pa  Po  d 

where  ^  =  linear  absorption  coefficient,  ^  is  the  linear  scattering 
coefficient  MT  =  Ms  +  «  total  attenuation  coefficient,  M  is  the  mass 

of  atmosphere  in  a  vertical  column  of  unit  cross  section  and  d  cm  in 

length,  p  is  the  atmospheric  density  at  altitude  h  and  p0  is  the  density 
at  N.T.P.  conditions. 

The  total  flux  at  altitude  h  produced  by  bremsstrahlung  photons 
of  energy  Ep  emitted  from  a  source  at  altitude  h  -t  d  is 

r  »  R 

=  /  dl 


The  ionization  produced  at  h  by  the  absorption  of  photons  with  energ 


E  is 
P 

£.1000 

p  35 
f  0 

The  total  ionization  rate  at  an  altitude  h  is 


tiy 


q(h'  V  ‘  '  Ep  \  Mr>  sec 

0 


qh  =  2  Ep) 
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where  q  is  summed  over  the  photon  energy  spectrum  for  the  bremsstrahlung 
produced  by  primary  electrons,  in  a  particular  energy  interval  Ee , 
absorbed  at  an  altitude  h  +  d;  and  then  summed  over  the  entire  spectrum 
of  primary  electrons,  absorbed  over  the  entire  altitude  range  of  interest. 

The  results  are  presented  in  Figure  III-15. 


FIG.  111-15  PEAK  BREMSSTRAHLUNG  IONIZATION  RATE  PROFILE 


The  combined  electron-bremsstrahlung  ionization  rate 
peak  primary  electron  flux  is  presented  in  Figure  1 1 1 —1 6 . 


profile  for 
The  temporal 
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FIG.  111-16  PEAK  IONIZATION  RATE  PROFILE  -  ELECTRONS  PLUS  BREMSSTRAHLUNG 

variation  of  the  ionization  rate  q  presented  in  Figure  1 1 1 -17  is 
a  smoothed  version  of  the  bremsstrahl ung  curve  recorded  by  Brown  and 
Campbell  (1962). 


4 .  Electron  and  Ion  Densities 


Knowing  the  ionization  rate,  q,  as  a  function  of  time  and  altitude, 

the  free  electron  concentration,  N  ,  and  the  ratio  of  negative  ions  to 

e 

electrons,  X,  can  be  computed  by  utilizing  the  following  well-known 
continuity  equations  (Mitra,  1952;  Whitten  and  Poppofi ,  1961) 
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FIG.  111-17  TEMPORAL  VARIATION  OF  IONIZING  PULSE 


where  is  the  dissociative  ion-electron  recombination  coefficient,  Qr^ 
is  the  ion-ion  mutual  neutralization  coefficient,  b  is  the  electron- 
oxygen  attachment  coefficient,  pis  the  negative-ion  photodetachment 
rate,  y  is  the  negative-ion  collisional  detachment  coefficient,  m  is 
the  concentration  of  the  attaching  species  (02  or  0),  and  n  is  the 
concentration  of  the  species  causing  col li sional-ty pe  detachment. 

Unfortunately,  the  computations  of  electron  density  profiles  could 
not  be  completed  in  time  for  this  report.  Calculations  of  temporal 
variations  in  electron  concentration  were  made  for  the  altitude  range 
75  to  120  km  and  are  presented  in  Figures  II 1-18  and  III-19.  Calculation 
of  an  electron  density  profile  was  made  for  the  peak  ionization  rate 
(51.6  minutes)  for  the  range  40  to  120  km  (Figure  III —20 ) ;  this  can  be 
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Altitude 


compared  with  the  peak  ionization  rate  profile  (Figure  III-16).  The 
coefficients  for  recombination,  neutralization,  attachment  and  detachment 
were  taken  from  the  appropriate  sections  of  this  report  The  initial 
conditions  for  X  were  also  taken  from  this  report  (see  section  III  A). 

The  coefficients  used  for  altitudes  below  80  km  become  increasingly 
more  speculative  as  the  altitude  decreases.  The  assumptions  made  in 
selecting  coefficients  for  the  lower  altitude  regions  are  noted  in 
Table  1 1 1-2. 

Although  the  observations  by  Brown  and  Campbell  (1962)  were  made 
in  the  vicinity  of  local  midnight,  calculations  indicate  that  the  atmos¬ 
phere  was  illuminated  throughout  the  observing  period  above  at  least 
3  km  Therefore,  there  should  be  no  question  about  the  occurrence  of 
photodetachment  above  the  ozone  layer. 

In  order  to  facilitate  the  calculations  and  provide  rough  estimates 
of  electron  concentrations  for  this  report,  simplifying  assumptions  were 
made  in  the  continuity  equations.  For  computations  of  electron  concentra¬ 
tions  above  75  km  ,  X  wns  considered  constant  with  respect  to  time.  For 
computations  below  75  km  ,  and  were  considered  equal .  The  results 
for  altitudes  above  75  km  are  considered  to  be  good  estimates  in  terms 
of  the  present  state  of  knowledge.  The  results  for  altitudes  below  75  km 
are  much  less  defensible  at  this  time  and  should  be  considered  to  be 
very  provisional.  The  assumptions  will  be  more  carefully  evaluated  and 
better  results  should  be  obtained  in  the  near  future  Time-dependent 
solutions  for  X  will  also  be  attempted. 

Comparisons  of  the  temporal-variation  curves  of  electron  densi¬ 
ties  (Figure  III-18)  and  the  temporal-variation  of  cosmic  noise  absorp¬ 
tion  (bottom  of  Figure  III-18)  are  noteworthy.  Inasmuch  as  cosmic  noise 
absorption  is  dependent  on  free  electron  densities,  the  siiilarity 
between  the  shape  of  the  electron  concentration  curves  and  the  shape  of 
the  cosmic  noise  curve  indicates  that  the  auroral  ionization  model  is 
at  least  partially  valid.  It  is  interesting  to  note  that  the  decay  of 
the  electron  density  curves  (after  the  two  ionization  maxima  at  51  6 
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Table  III-2 


INITIAL  CONDITIONS  AND  RECOMBINATION  PARAMETERS 


Recombination  Parameters 

Initial 

Condi tions 

(km) 

aD 

a 

<*i 

b 

pm 

c 

d 

p 

r 

r  n 

e 

N 

e 

l 

cm 

/ 

sec 

cm 

sec 

sec 

-1 

sec” 1 

sec 

-  i 

cm^ 

3 

A 

25 

2.0 

X 

10"6 

2.1 

X 

10"7 

3.6 

X 

10  4 

0.1 

1  6 

X 

io-2 

<i 

3  6 

x  10  s 

30 

2  0 

X 

10"6 

9  3 

X 

10"8 

8  1 

X 

103 

0.1 

1.5 

X 

10"2 

<i 

8  0 

x  10 4 

35 

2.0 

X 

10"G 

4  5 

X 

10"8 

1.7 

X 

10  3 

0  1 

8.5 

X 

10"  3 

i 

2  0 

x  104 

40 

2  0 

X 

10 -G 

1  8 

X 

IO-8 

4  1 

X 

102 

0  1 

6.0 

X 

IO"3 

8 

4 

x  103 

45 

2  .0 

X 

10"G 

1  0 

X 

10  “R 

9  7 

X 

10  s 

0  1 

4.0 

X 

10“  3 

20 

9  7 

x  102 

50 

2.0 

X 

10-G 

1  0 

X 

10“8 

2.5 

X 

101 

0.1 

6.0 

X 

10"  3 

30 

250 

55 

2.0 

X 

10“G 

1  .0 

X 

10“8 

9.5 

X 

10° 

0.1 

4  0 

X 

10-3 

50 

95 

60 

2.0 

X 

10"G 

1  0 

X 

IO"8 

2  4 

X 

10° 

0  1 

9  0 

X 

10-4 

2  .x 

102 

24 

65 

1 .5 

X 

10"6 

1  0 

X 

10“8 

6  9 

X 

10-1 

0  1 

4  0 

X 

10-4 

2  5  x 

102 

6  9 

70 

1  0 

X 

10"® 

1  0 

X 

10"* 

2  .0 

X 

10"1 

0  1 

2  0 

X 

IO-3 

4  x 

102 

2  0 

75 

8  0 

X 

10"7 

1  .0 

X 

10“8 

4  5 

X 

10"2 

0  1 

4  .0 

X 

10-3 

6  x 

102 

0  3 

80 

6  0 

X 

10~7 

1  .0 

X 

10“8 

9  0 

X 

10-3 

0  1 

1  0 

X 

10-2 

6  x 

102 

4.5 

x  10"2 

85 

5.1 

X 

10-7 

1 .0 

X 

10~8 

2.0 

X 

10"3 

1  0 

1.5 

X 

10-2 

1 .5  x 

10  3 

2 

x  10“3 

90 

5  0 

X 

10"7 

1  .0 

X 

10'  8 

4.0 

X 

10"  3 

1.0 

3.0 

X 

10“2 

2  x 

10  3 

4 

x  IO-3 

95 

4  .7 

X 

IO1"7 

1  .0 

X 

10“* 

6.7 

X 

10-3 

1.0 

5.0 

X 

10-2 

6  x 

10  3 

7 

x  IO-3 

100 

4.2 

X 

10~7 

1  .0 

X 

10"8 

2.7 

X 

io-3 

1,0 

2.0 

X 

10~2 

3  x 

104 

2.7 

x  10-3 

105 

3.8 

X 

10-’ 

1  0 

X 

IO'8 

2.2 

X 

io-3 

1  0 

1.6 

X 

10-2 

10: 

2  2 

x  IO"  3 

110 

3.4 

X 

10"7 

1  .0 

X 

10'8 

1  .3 

X 

10"3 

1.0 

1  .0 

X 

10-2 

2  x 

10s 

1.2 

x  IO"  3 

115 

3  0 

X 

10-’ 

1.0 

X 

IO-8 

5.4 

X 

10-4 

1  .0 

4  .0 

X 

IO"3 

2  x 

10s 

5.4 

x  10~4 

120 

2  6 

X 

IO"7 

1.0 

X 

10"8 

2.7 

X 

10”4 

1  0 

2.0 

X 

10' 3 

2  x 

10s 

2.7 

x  10"4 

Values  of  ’I}  below  75  km  assume  a  clustering  of  molecules,  values  of  a0  above 
80  km  are  from  curves  of  temperature  dependence  in  this  report,  Chapt  IV 

Three-body  Thompson  coefficient  below  45  km. 

Above  K5  km  radiative  attachment  to  0  considered  dominant,  below  80  km  three- 
body  attachment  to  02  considered  dominant 

Above  85  km  detachment  from  o"  considered  dominant,  below  80  km  detachment 
from  02  considered  dominant 

Above  05  km  associative  detachment  considered  dominant,  below  GO  km 
coll isional  detachment  considered  dominant. 

Above  85  km  initial  )  considered  equal  to  .0  0~ 

Below  80  km  initial  X  considered  equal  to  ( 02 |  02“ 

where 


hr 

[0] 

fs, 

02  , 
p(0‘) 

i-  (,02  ) 


radiative  attachment  coefficient  for  0 
concentration  of  atomic  oxygen 

three-body  attachment  coefficient  for  02  =  K[02j  (see  Sec.  Ill  A  3 
concentra ti on  of  molecular  oxygen 
photo  detachment  rate  for  o” 
photo  detachment  rate  for  02~ 


and  75  mins)  reassemble  the  decay  of  the  cosmic  noise  absorption  curve 
more  closely  as  the  altitude  decreases.  The  results  are  not  yet 
sufficiently  complete  to  determine  the  altitude  range  that  is  most 
effective  in  producing  this  particular  cosmic  noise  record.  However, 
the  data  do  seem  to  indicate  that  the  lower  altitude  electron  densities 
are  important.  Inasmuch  as  the  collision  frequency  increases  rapidly  with 
decreasing  altitude  and  the  electron  density  in  the  range  75  km  to  60  km 
decreases  slowly  with  altitude,  (and  in  fact  increases  to  a  maximum  point 
at  60  km  according  to  Figure  III-20)  it  seems  reasonable  to  expect  that 
the  electrons  in  this  region  are  very  important  in  producing  cosmic  nosie 
absorption . 

Ionization  in  this  region  is  caused  by  bremsstrahlung  rather  than 
primary  electrons.  (Compare  Figures  III-13  and  III-16,  Therefoie, 
if  the  apparently  important  influence  of  lower  altitude  electron 
concentrations  is  confirmed  by  future  results,  the  importance  of 
bremsstrahlung  will  also  be  demonstrated. 

D .  Chemionization 

Chemionization,  which  is  the  production  of  electron-ion  pairs  by 
means  of  chemical  reactions,  has  been  investigated  in  the  laboratory 
only  relatively  recently  (Kunkel  and  Gardner,  1955;  Kunkel ,  1956;  Smith 
and  Gatz,  1961).  Since  various  reactants  which  can  produce  chemioniza¬ 
tion  are  present  in  the  upper  atmosphere,  it  is  considered  worthwhile 
to  investigate  the  possible  effects  on  ionization  levels  and  relaxation 
times.  Several  energetically  possible  processes  are: 


N  +  N  +  NO  -*  N0+  +  N2  +  e  (l) 

N  +  N  +  M  ->  Nz+M  (2a) 
Nj  +  NO  -<  N2  +  NO  +  e  i2b) 
Na  +  0  +  0  -  Na+  +  02  +  e  (3, 


where  M  represents  an  unspecified  third  body.  Processes  (l)  and 
are  exothermic  and  (3)  is  endothermic  by  an  amount  approximately  equal 
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to  the  mean  kinetic  energy  of  the  reactants.  In  process  (2)  the  inter¬ 
mediate  excited  state  N2  can  also  decay  by  emission  of  radiation.  Hence 

it  must  be  metastable  against  radiation  if  the  collisional  mode  (2b) 
domi nates. 

The  corresponding  electron  production  rates  q  are 


q,  =  kj  S’  j 2 

NO) 

(■0 

„  k2a  kofc 

N  2  ( NO ] [ M ] 

q2  -  j— 

(5) 

~  + 

T 

k2b  1K0J 

da  =  k a  i Na j 

Oj2 

(6) 

where  the  ka’s  are  the  rate  constant-;  of  the  corresponding  processes 
listed  above  and  r  is  the  lifetime  of  N*  with  respect  to  radiation. 

The  concentrations  of  the  indicated  atomic  and  molecular  species  at  an 
altitude  of  80  km  are  estimated  to  be  of  the  order 


10*  cm-3  (Barth,  1981) 


N0j 

106  cm-3 

(Barth,  19(11) 

Oi 

It)12  cm-3 

(Barth,  1961) 

M 

10* 4  cm” 3 

(Minzner,  et  al . ,  1959; 

Na 

104  cm-3 

(Bates,  19(10) 

Actually,  Bates  reports  several  measurements  of  the  number  ol  sodium 
atoms  m  a  vertical  column  tc  be  of  the  order  2-5  x  109  cm"2.  Since 

most  of  the  sodium  lies  in  the  altitude  range  80  to  100  km  with  n  peak 

at  85  km,  the  concent ra t ■ on  of  sodium  at  85  km  is  estimated  to  be 
104  cm"3  if  we  assume  that  all  of  it  exists  as  free  atoms. 

Assuming  (l  that  the  reaction  rates  are  of  the  order  JO"10  cm3  sec"1 

for  two-body  processes  and  lO’30  cm6  sec’ ‘  for  three-body  processes 
(these  values  are  upper  limits  for  the  rate  constants),  (2)  that  the 
radiative  lifetime  of  N2*  is  1  sec,  and  (3)  that  no  activation  energy 
is  required  in  the  reactions,  we  obtain 
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q,  ~  10“ 12  cm-3  sec-1 
q2  ~  10”8  cm-3  sec-1 
q 3  ~  10“2  cm-3  sec-1 


of  which  q,  and  q2  are  so  small  that  equilibrium  concentrations  of  elec¬ 
trons  produced  thereby  are  only  of  the  order  10  2  and  1  cm  ,  respectively; 
we,  therefore,  neglect  them.  However, process  (3)  may  be  of  marginal 
significance  as  a  means  of  producing  electrons  at  night  in  the  region 
bO  to  90  km.  In  order  to  determine  the  nighttime  contribution  to  the 
electron  concentration  from  chemioni zation  one  must  solve  a  number  of 
first-order  nonlinear  differential  equations  with  and  without  the  chem- 
ionization  source  term;  the  equations  containing  the  source  term 


k3  rNa  ,0  are 

^  =  k3  Naj  0  j 2  -  aD(N0+)[K0+][e]  +  7  [0] [0^]  -  fr(0t ) '02 ] [e  ]  -  K[02]2[e] 

(?) 

d  Xa- J  =  k  3 i Na  0]2  -  Q  (Na+ ,  02)[Na+][02J  (8) 

dt  i 

‘‘  =  -  ai(N0+,  02  )[.NO+][ojj  -  aD(N0+);N0+  .e  (9) 


=  K  02  j2  ,e  +  P(03  )  [02  .  ,e  j  -  I  Q.  (Na  +  ,  02  )  [Na+  +  Cl^  (NO  ,  02  NO  I)  [02  ] 

-  y'Ojlojj  (10) 

where  a  and  a  are  the  dissociative  and  ion-ion  recombination  coefficients 
D  l 

of  the  species  indicated,  K  is  the  three-body  and  (j  the  radiative  electron 
attachment  coefficient  of  02 ,  and  y  is  the  associative  detachment  coeffi¬ 
cient  of  02 .  Actually  negative  ions  other  than  02  and  positive  ions  other 
than  N0+  and  Na+  may  also  bo  important.  Their  inclusion  would  further 
complicate  an  already  formidable  set  of  equations.  Not  only  is  the 
solution  difficult  but,  more  important,  the  values  of  many  of  the  rate 
constants  are  unknown. 
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One  can  arrive  at  rough  estimates  of  the  electron  concentration  at 
several  hours  after  sunset  by  adding  equations  (7)  and  (10)  and  assuming 

that  the  negative  ion-electron  concentration  ratio  X  =  is  nearly 

constant  in  time.  G- 

d[e]  q 3  ai  (Na+,02 ) [Na+  +  (NO+,oj) [NO+ ] 

dt  =  77\  "  *  - - — — - - - — - 


We  further  assume  that 


qd(no+)[no+J 

[Na+  +  [NO+] 


[e]2  (11) 


X  *=  1 

(NO  ,  02 )  ^  10  8  ^this  paper) 
a.  (Na+,  02)  =  10"*  (Bates,  1960) 
ap(NO  )  =  10  7  (this  paper) 

initial  *  1300 

Electron  concentrations  at  various  times  after  sunset  were  computed 
using  the  values  of  the  coefficients  and  concentrations  given  above 
for  two  cases;  (1)  without  the  occurrence  of  chemi oni zation,  and  (2) 
with  the  occurrence  of  chemi oni zation.  The  results  are  presented  in 
Table  IV-6 . 


Table  IV-6 

NIGHTTIME  ~  .ECTRON  CONCENTRATIONS 
(at  85  km) 


Time  after 
Sunset 
(hrs ) 

e.  Without 
Chemioni zation 
(cm-3) 

.ej  With 

Chemi oni zation 
(cm'5) 

a 

effective* 
(cm3  sec"1 ) 

1 

900 

1000 

6.5  x  10-8 

2 

670 

800 

6.4  x  10-8 

3 

540 

690 

5.5  x  10-8 

4 

440 

600 

5.5  x  10"8 

5 

380 

540 

5.0  x  10"8 

6 

340 

500 

4.2  x  10-8 

7 

300 

470 

3  x  10-8 

♦With  chemionization  occurring  but  not  included  in  the 
calculation  of  a 

effective’ 
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It  is  evident  that  even  if  the  rather  extreme  value  assigned  to 
the  electron  production  rate  q3  is  correct,  the  effects  of  chemioniza- 
tion  on  the  nighttime  ionosphere  are  marginal  at  most.  If  chemioniza- 
tion  is  present  and  is  not  considered  in  deriving  the  effective  night¬ 
time  recombination  coefficient  (from  observation),  the  latter  will 
appear  to  decrease  steadily  throughout  the  night;  this  is  shown  in  the 
right-hand  column  of  Table  IV-6.  Mitra  (1957)  has  reported  such  an 

observed  decrease  in  a  between  85  and  110  km  altitude.  He 

effective 

has  attributed  it  to  electron  attachment.  However,  the  lifetime  of 
electrons  with  respect  to  attachment  is  only  of  the  order  of  a  few 
minutes  in  this  region  and  the  attachment  process  thus  approaches 
equilibrium  rapidly.  Nevertheless,  one  cannot  consider  this  evidence 
as  definitely  supporting  the  suggestion  of  chemionization ,  particularly 
since  process  (3)  cannot  be  expected  to  be  effective  above  95  km.  The 
question  is  open. 

Chemionization  is  of  potential  importance  in  the  ionosphere  only 
»hen  other  sources  of  ionization  are  absent  or  very  small.  For  example, 
the  ultraviolet  radiation  from  a  nuclear  explosion  is  expected  to  cause 
momentarily  a  very  high  degree  of  dissociation  of  Nz  and  02  in  the 
vicinity  of  the  fireball,  perhaps  10/C.  Large  concentrations  of  NO  will 
also  form  and  the  electron  production  rate  from  chemionization  will  be 
quite  large  a  value  of  q,  =  106  cm-3  sec-1  would  not  be  unreasonable. 
However, the  local  photoionization  will  be  so  large-as  to  completely  mask 
the  chemionization  effect.  After  the  initial  flash  of  ultraviolet 
radiation  and  X-rays,  the  atomic  nitrogen  and  oxygen  concentrations 
will  rapidly  decay  to  their  normal  values  since  radiation  from  fission 
fragments  is  not  very  effective  in  causing  dissociation.  Therefore  one 
would  expect  chemionization  to  cease  immediately  after  the  initial 
radiation . 
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IV  DE-IONIZATION  PROCESSES 


A. 


Negative  ions 


Li  uuuc  t  ion 


knowledge  of  .he  nature  of  the  nee, ml  species 

*"  th°  the  possible  specie,  of 

negative  ion,  o'.  <.  oj.  n, .  ,„d  »hlc„  ^ 

negative  is  a  question  which  ean  only  be  definitely  answered 

by  mass  spectrometric  measurement®  It  mau 

urement®.  It  may  even  turn  out  that  some 

species  not  considered  heretofore  Is  the  most  important  one.  Having  no 

idea  of  the  nature  of  such  a  species,  however,  we  must  of  necessity 

-estrict  out  discussion  to  the  six  previously  mentioned. 

..  „  ,ThC  nUmCr°US  Pr°CCSSeS  inVOlvi"«  formation,  destruction,  and 
shuffUng  or  interchange  of  negative  ion  species  may  be  broadly  classi- 
fied  as  three-body  attachment,  radiative  attachment,  collisional  type 
de.aeh„„„t,  Pho, edetachaten,  ,  i«.  reco.b, nation ,  cha^c  transfer. 

“  T"  *°  d“s"  The  <‘“s  npt  —  ■“  «*- 

Possible  b„,  do.,  ,„ei„de  processes  CPns,derea  „„ 

(except  diffusion  to  dust): 


Three-body  attachment 


°2  +  02  +  e  -  oj  ♦  02 
°2  +  0  +  e  —  o2  +  0 
02  +  0  +  e  02  +  o” 

Radiative  attachment 

°2  +  e  -  Oj  +  hv 

0  +  e  o'  +  hv 

H  +  e  -  H~  +  hv 


(la) 

(lb) 

(lc) 


(ld) 

(le) 

(lf) 
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Photodetachment 


02 

+  hv 

—  02  4  e 

(lg) 

o' 

+  hv 

—  0  4  e 

(lh) 

h" 

+  hv 

-  H  4  e 

(li) 

NOj 

,  +  h\ 

—  N02  4  e 

(lj) 

oh" 

+  h\ 

-  OH  4  e 

(lk) 

°3 

+  hv 

-  03  4  e 

(10 

Col li sional-ty pe  detachment 

Cl 

4  M  — 

02  4  M  4  e 

(lm) 

0: 

4-  0  • 

03  ♦  e 

(In) 

02 

4  0- 

O3  +  hv;  O3  +  0  —♦  Oj  +  O2  +  0 

(lo) 

O" 

4  0- 

02  4  e 

(IP) 

if 

4  H  - 

Hj  4  e 

(iq) 

OH_ 

4  H 

-<  H20  4  e 

(lr) 

N  + 

no2 

■4  N2  4  02  4  e 

(Is) 

0* 

4  0  - 

0  2  4  0  2  4  e 

(It) 

Ion-Ion  recombination  of  the  types 

02  ,  02  — »  02  +  02  (I*1) 

0"  +  N2  NO  -t  N  (iv) 

02+02+M-*02+02+M  ( lw ) 

Charge  transfer  and  ion-atom  interchange 

02  +  0  -<  02  +  o" 
o"  +  NO,  -  0  +  N02 
02  +  H  — »  02  +  H 


(lx) 

(iy) 

(lz) 
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+  OH  02  +  OH 

(laa) 

+  OH  0  +  Otf 

(lbb) 

+  N2  -  N  +  N0~ 

(lcc) 

■  02  +  M  -  NO 2  +  M 

(ldd) 

+  02  +  02  -  02  +  02 

(lee) 

+  0  +  02  -  Oj  +  02 

(Iff) 

2.  Electron  Attachment 


The  formation  of  02  by  the  three-body  process  (la)  has  been  investl 
Rated  by  Chanin  and  coworkers  (1959)  who  found  the  rate  constant  to  be 
sufficiently  larKe  to  make  the  reaction  of  preat  potential  importance 
in  the  D-reRion.  For  example  at  a  temperature  of  -  200°K  which  is 
typical  of  D-reRion  altitudes  one  obtains  a  value  for  the  rate  constant 
of  K  _  (1.4  t  0.5)  x  10  30  cm3  sec.  These  investipators  also  found  that 
N2  Is  not  an  efficient  "third  body"  in  this  type  of  reaction.  beinR 
only  about  2i  as  effective  ns  0,.  The  effectiveness  of  atomic  oxypen 
as  a  "third  body"  has  not  been  i nvestiRated.  However,  ills  reasonable 
to  assume  that  it  is  not  more  efficient  than  02  and  we  therefore  noRlect 
process  (lb)  because  of  the  small  (relative  to  02  )  concentrations  of 
atomic  oxypen  in  the  D-laycr. 


OwinR  in  part  to  obvious  experimental  difficulties,  three-body 
attachment  to  atomic  oxyRen,  process  (1c),  has  never  been  investiRated. 
However ,  we  may  arrive  at  a  very  rouRh  estimate  of  the  ratio  of  the 
rate  constant  to  that  of  process  (lb)  from  the  work  of  Curran  (1961) 
on  dissociative  attachment  to  ozone.  Curran  found  that  the  branchinR 
ratio  of  the  processes 


+  0 

(2a) 

+  O2 

(2b) 

is  about  7/11.  If  we  assume  that  the  rate  of  electron  attachment  to 
02  and  0  is  about  the  same  in  the  unbound  state  (processes  (lb)  and 
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(lc))  as  in  the  bound  state  (process  (2)),  we  can  assign  a  branching 
ratio  to  processes  (lb)  and  (lc)  of  order  unity.  We  have  already 
advanced  arguments  for  neglecting  process  (lb)  and  thus  we  neglect 
(lc)  also. 

Branscomb  and  coworkers  at  the  National  Bureau  of  Standards  (Branscomb, 

Burch,  Smith  and  Geltrr.an,  1958;  Burch,  Smith  and  Branscomb,  1958;  Smith 

and  Burch,  1959;  Smith  and  Branstomb,  1960)  have  measured  the  photo- 

detachment  cross  sections  of  0  ,  02 ,  and  H  for  photon  energies  near 

threshold  up  to  several  electron  volts  (to  be  discussed  in  detail  later). 

The  radiative  attachment  cross  section  a  (0  )  can  be  easily  obtained 

3tt 

from  the  photodetachment  data  by  use  of  the  principle  of  detailed 
balancing : 

-  £>■  r,  <3) 

where  <  is  the  momentum  of  the  incident  photon  divided  by  if,  k  is  the 

K_ 

momentum  of  the  incident  electron  divided  by  -K,  —  is  the  ratio  of  the 

K° 

statistical  weight  of  the  0  ground  state  to  that  of  the  oxygen  atom 
ground  state,  and  <Jdct  is  the  photodetachment  cross  section.  One  easily 
obtains  the  effective  radiative  attachment  rate  constant  for  any  electron 
temperature  by  averaging  the  product  of  attachment  cross  section  and 
electron  speed  over  all  electron  energies: 


e 


r 


o  E  e_E  kT  dE  (1) 
att  ' 


where  T  is  the  electron  temperature,  m  the  electron  mass,  E  the  electron 

energy,  and  k  is  Boltzmann’s  constant.  Branscomb  and  coworkers  (1958) 

found  that  the  radiative  attachment  coefficient  for  atomic  oxygen  is  a 

slowly  varying  function  of  electron  temperature,  decreasing  from 

£  =  1.3-i  x  10" 15  cm3  sec  at  an  electron  temperature  of  250°  to 

r 

£  =  1.25  x  10-is  cm3/ sec  at  a  temperature  of  1000  K.  Actually  use  of 

the  principle  of  detailed  balancing  for  computation  of  a  is  not 
fully  justified  for  0  in  view  of  recent  evidence  (Schulz  ,  1962)  that 

lower-lying  energy  states  of  the  0  ion  may  also  exist.  Hence  the  value 
of  P  given  here  must  be  recognized  as  provisional  only. 
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In  their  experiment  on  02  photodetachment  cross  sections  Burch 
2.\  ‘Li." »  (1958)  found  that  the  state  of  the  02  ion  after  attachment  has 
gerade  symmetry,  thus  permitting  one  to  express  the  energy  dependence 
(Geltman ,  1958)  of  the  cross  sections  as 


det^°2^  =  E(E  ”  E0/3  2  A0  +  Aj  (E  -  E0)  +  A2(E  -  E0)2  +  ••■] 

(5) 

where  E  is  photon  energy,  E0  is  threshold  energy,  and  the  A  are  constants 

n 

related  to  the  electric  dipole  transition  matrix  elements.  In  practice 
the  A^  are  uncalculable  to  any  acceptable  degree  of  precision  and  must 
be  found  empirically.  The  first  two  were  found  to  be  A0  =  0.370  x 
10-18  cm2  eV-*/2,  A,  =  -0.071  x  10"18  cm2  eV7'2.  The  0~2  created  In 
the  experiment  apparently  underwent  very  few  collisions  (of  the  order 
of  one  to  ten  per  ion)  (Smith  and  Branscomb,  1960)  prior  to  electron 
detachment.  We  shall  cite  evidence  later  which  seems  to  indicate  that 
the  02  is  formed  in  an  excited  state  and  then  deactivated  only  after 
many  collisions  (»  10).  Assuming  that  this  interpretation  is  correct, 
the  02  ions  must  have  remained  essentially  unchanged  during  the  course 
of  the  experiment. 


Invoking  the  principle  of  detailed  balancing,  one  obtains  the 
following  equation  for  o  with  the  aid  of  (5) 


att 


*  (£  +  E°) 


3  12 


2mc4 


An  +  Alt.  +  A2t;  2  + 


■  J  —  (fi) 

Ko 


whore  c  is  the  kinetic  and  me2  the  rest  energy  of  the  electron.  Equa¬ 
tions  (4)  and  (6)  together  with  the  previously  stated  values  of  A0  and 

A,  yield  a  value  for  ^(0^)  of  ~  10"20  cm»/scc  at  an  electron  tempera¬ 
ture  of  250°K . 


Of  course,  the  principle  of  detailed  balancing  is  not  really 
applicable  here  because  neither  02  nor  02  is  necessarily  left  in  the 
ground  state  after  electron  attachment  and  detachment;  excited  vibra¬ 
tional,  rotational,  and  possibly  electronic  states  can  be  occupied.  In 
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other  words,  the  process  is  not  reversible  and 


l<  O2,  hv  |  H1  |  02,  ej  f  t 


02 ,  e  H1  |  02,  hV  >  |  , 


(7) 


where  H1  is  the  electric  dipole  operator.  However,  it  is  thought 

unlikely  that  the  effect  of  such  allowed  excited  states  would  increase 

the  value  of  p^  by  more  than  two  to  three  orders  of  magnitude  from 

that  given  above.  In  order  to  be  of  any  importance  in  determining  the 

rates  of  ionospheric  processes  we  must  have  p  >  10“17  cm3/sec.  Hence 

r  ' 

it  is  unlikely  (though  not  certain)  that  radiative  attachment  of  elec¬ 
trons  to  02  is  an  important  atmospheric  reaction  under  normal  conditions; 
this  process  may  be  of  marginal  significance  at  altitudes  of  85  km  and 
above.  It  can  be  shown  that  the  radiative  attachment  coefficient 
increases  with  increasing  temperature.  At  an  electron  temperature  of 
1000°K ,  for  example,  P^)  >  10“19  cm3/ sec,  an  increase  of  about  a 
factor  of  ten  over  the  value  at  a  temperature  of  250°K.  Hence  under 
disturbed  conditions  (e.g.,  SID  or  nuclear  blackout)  radiative  attach¬ 
ment  to  02  may  conceivably  be  of  some  importance.  Further  experimental 
investigation  is  definitely  in  order. 


Smith  and  Burch  (1959)  have  measured  the  photodetachment  cross 
section  of  H  for  various  electron  energies  above  threshold.  Use  of 
the  principle  of  detailed  balancing  (which  can  be  justified  as  in  the 
case  of  0  yields  a  radiative  attachment  coefficient  of  -v  10",G  cm3  sec 
at  a  temperature  of  200 °K.  Because  of  the  relatively  small  concentra¬ 
tions  (n(H)  ^  109  cm"3)  (Handbook  of  Geophysics  ,1960 )  direct  attachment 
to  this  species  proceeds  at  a  much  smaller  rate  than  direct  attachment 
to  atomic  and  molecular  oxygen  and  hence  process  (If)  can  be  neglected. 


We  are  thus  left  with  processes  (la)  and  (le)  as  the  dominant  direct 
attachment  reactions.  At  an  altitude  of  85  km  the  rate  of  electron  attach¬ 
ment  to  both  atomic  and  molecular  oxygen  is  about  ~  10-3  sec-1 ,  assuming 
concentrations  of  [oa]  O  x  10”  cm*3  and  [0]  ~  5  x  10“  cm'3  (Hand¬ 
book  of  Geophysics,  1960).  Below  this  altitude  the  rate  of  the  three-body 
mode  (la)  increases  rapidly  due  to  increasing  concentration  of  02 ,  while 
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the  rate  of  process  (le)  decreases  slowly.  Hence  we  conclude  that  radi¬ 
ative  attachment  of  electrons  to  atomic  oxygen  is  important  only  at  the 
top  of  the  D-region  and  that  process  (la)  provides  the  dominant  attach¬ 
ment  mechanism  below. 

It  should  be  noted  that  we  have  excluded  direct  attachment  to  N02 . 
Curran  (1962)  has  apparently  established  this  process  to  be  forbidden. 

3 .  Electron  Detachment 

Some  of  the  electron  detachment  processes  (lg,  lh,  li,  lm)  have 
also  been  investigated  recently  but  the  results  have  not  been  sufficient 
for  the  development  of  a  reliable  model  of  electron  detachment  in  the 
ionosphere . 

Photodetachment  cross  sections  in  the  photon  energy  range  0  to 
3  eV  have  been  measured  for  o"  and  02  by  Branscomb  and  coworkers 
(Branscomb,  Burch,  Smith  and  Geltman,  1958;  Burch,  Smith,  and  Branscomb, 
1958)  at  the  National  Bureau  of  Standards.  The  o“  cross  section  was 
found  to  be  a  rapidly  increasing  function  of  photon  energy  between 
threshold  (EA  =  1.465  eV)  and  ~  1.6  eV  (a  =  5  x  10"1®  cm2)  after  which 
it  levels  off,  reaching  7  x  10-1®  cm2  at  a  photon  energy  of  ~  3  eV. 

These  results  are  in  good  agreement  with  the  theory  of  photodetachment 
from  negative  atomic  ions  developed  by  Klein  and  Brueckner  (1958).  The 
ionospheric  specific  photodetachment  rate  is  obtnined  from  the  equation 

P  =/  CTdet<f>  dX  (8) 

where  is  the  photodetachment  cross  section  and  I(X)  is  the  solar 

photon  number  flux  at  wavelength  incident  on  the  ionosphere.  The 
result  for  0  is  p  =  1.4  sec"1  (Smith,  Burch,  and  Branscomb,  1958). 

Because  of  the  possible  presence  of  excited  states,  the  situation 
for  photodetachment  from  02  is  more  complicated;  one  would  in  general 
expect  to  find  different  cross  sections  for  detachment  from  ground  and 
excited  states,  partly  because  of  the  difference  in  predicted  threshold 
energies  and  partly  because  of  different  transition  matrix  elements 
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for  the  two  cases.  Geltman  (1958)  has  computed  the  energy  dependence 
near  threshold  of  the  photodetachment  cross  section  of  0~2  and  has 
obtained  relation  (5).  In  their  investigation  of  the  problem,  Burch 
and  coworkers  (1958)  were  able  to  determine  A0  and  A,  as  well  as  E0 
by  a  least  squares  fit  of  their  experimental  data  to  the  theoretical 
form  of  det.  The  value  of  the  threshold  energy  (E0  =  0.15  eV)  was 
considerably  smaller  than  that  obtained  by  Phelps  and  Pack  (1961) 

(E°  =  0,46  eV)  by  m°anS  °f  Visional  detachment  experiments.  This 
evidence  seems  to  indicate  that  <£  is  formed  in  an  excited  (probably 
vibrational)  state.  The  gas  pressures  and  ion  residence  times  in  the 
swarm  experiments  of  Phelps  and  Pack  were  sufficiently  large  to  permit 
colllsional  de-excitation  to  occur;  they  estimated  that  approximately 
10  collisions  of  each  02  with  neutral  molecules  occurred  during  the 
experiment.  On  the  other  hand,  the  gas  pressures  as  well  as  the  resi¬ 
dence  times  in  the  photodetachment  experiments  (Smith  and  Branscomb, 

I960;  were  so  small  that  the  mean  number  of  collisions  of  each  0~2  ion 
with  a  neutral  particle  was  only  of  the  order  of  one  to  ten. 

In  contrast  to  the  results  of  the  experiments  mentioned  above 
chemical  experiments  (Pritchard,  1953)  seem  to  indicate  an  electron 
affinity  of  slightly  loss  than  1  eV.  A  spectrophotometric  measurement 
by  Jortner  and  Sokolov  (1961)  and  a  measurement  by  Curran  (1961)  using 
data  from  an  experiment  on  dissociative  attachment  to  ozone  yielded 
somewhat  smaller  values  for  the  electron  affinity,  ~  0.75  eV  and  -  0.58  eV, 
respectively.  It  is  difficult  to  reconcile  these  results  with  those 
of  Phelps  and  Pack  (1961).  A  possible  explanation  is  the  existence  of 
a  metastable  vibrational  state  at  E0  *  0.,6  eV,  the  de-excitation  from 
this  to  the  ground  state  being  of  a  high  order  of  forbiddenness.  m 
order  to  complete  this  explanation  one  must  assume  that  in  the  experi¬ 
ments  of  Pritchard,  Jortner,  and  Curran  <£  is  either  formed  in  the  low¬ 
ing  state  or  undergoes  relaxation  in  such  a  way  as  to  avoid  the  meta¬ 
stable  state. 
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It  Is  expected  that  most  of  the  0,  ions  in  the  sunlit  D-region 
will  be  in  the  ground  state  or  a  low-lying  metastable  excited  state 
because  of  the  relatively  high  gas  pressures  and  long  lifetimes  (>  2 
seconds)  of  the  ions.  If  one  approximates  the  solar  radiation  spectrum 
in  the  near  infrared  and  visible  regions  by  a  black-body  distribution 
and  substitutes  it  for  jjl  in  equation  (8),  one  obtains  an  approximate 
threshold  energy  dependence  of  p 

,  .  .  Igfl  -  E») 


£i=a 

HeH 


*=  C 


(9) 


For  a  solar  surface  temperature  of  6000°K  one  obtains 


(10) 


By  substituting  the  observed  solar  spectrum,  -^1  ,  in  equation  (8) 
Smith  £t  al.,  (1958)  obtained  a  value  of  0(0.15  eV)  =  0.44  sec-1. 
Making  the  reasonable  assumption  that  A0(0.15  eV)  >  Ao(0.46  eV),  we 
have  ^(0.46  eV,  0.2  sec" 1 .  Further  investigation  of  this  problem 
is  required  before  we  can  obtain  a  reliable  D-region  photodetachmont 
rate.  We  shall  provisionally  adopt  a  value  of  o(oj)  =  0.1  sec-1. 

Smith  and  Burch  (1959)  have  investigated  photodetachment  from  H_. 
Their  cross  section  data  together  with  the  solar  spectral  data  used 
for  the  computation  of  o(cf)  and  r(0~7)  (Smith,  Burch,  and  Branscomb, 
1958)  yield  a  detachment  rate  of  p  =  8  sec-'. 


Photodetachment  from  W2  ,  Off,  and  0-,  has  not  been  investigated 
although  some  idea  of  the  electron  affinities  is  available.  Curran 
(1962)  has  shown  that  the  electron  affinity  of  N0~  is  greater  than 
3.82  eV  by  causing  NO,  to  undergo  charge  transfer  with  Cl  (whose 
EA  =  3.82  eV).  The  ion  0H~  apparently  has  an  electron  affinity  lying 
somewhere  between  2  and  3  eV  (Field  and  Franklin,  1957)  while  that  of 
°3  is  estimated  to  be  of  the  order  of  3  eV  (Curran,  1961). 
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It  has  been  proposed  that  N02  may  be  the  dominant  negative  ion  in 
the  D-region.  This  suggestion  originated  with  the  apparent  observation 
by  Johnson  and  coworkers  (Johnson,  Heppner,  Holmes,  and  Meadows,  1958) 
of  negative  ions  of  mass  number  -46  in  the  E-region.  The  reliability  of 
these  measurements  is  uncertain;  the  Noj  may  have  been  a  product  of 
the  rocket  exhaust. 

It  is  known  that  just  after  sunset  and  just  prior  to  sunrise  when 
the  sun  s  rays  which  illuminate  the  D-region  must  first  pass  through 
the  ozone  layer  (Figure  IV-1 )  ,  the  electrons  in  the  D-region  rapidly 


FIG.  IV-1  SUNRISE-SUNSET  EFFECT 

At  sunset  the  neor  UV  solor  rodiotion  which  would  otherw.se  illuminote  thot  port 
of  the  D-region  lying  to  the  left  of  line  AB  is  obsorbed  by  the  ozone  loyer  when 
the  sun  drops  below  the  horizon  with  respect  to  on  observer  of  C.  Coincident 
with  the  obsorption  of  this  rodiotion  is  the  ropid  decrease  in  electron  density 
'Probably  due  to  electron  ottochment)  in  the  region  to  the  left  of  line  AB.  The 
opposite  effect  (ropid  detochment  of  electrons  when  the  sun’s  roys  poss  obove 
the  ozone  loyer)  occurs  ot  sunrise. 


attach  to  form  negative  ions  (sen,  for  example,  Reid,  1961).  This 
is  called  the  "sunrise-sunset  effect."  One  infers  from  this  that  the 
solar  radiation  responsible  for  photodetachment  must  be  strongly  absorbed 
by  ozone  and  must  therefore  lie  in  the  near  t'V  spectral  region.  If  one 
computes  the  spectral  detachment  rate  for  o"  and  h"  by  means  of  the 


equation 


where  the  symbols  have  the  same  significance  as  in  equation  (8),  one 
obtains  the  curves  shown  in  Figure  IV-2  .  The  near  UV  obviously  is  of 
little  importance  in  photodetachment  from  0~  and  H_.  We  conclude  from 
this  that  0  and  H  are  not  the  dominant  negative  ions  present  in  the 
D-region.  If  one  then  computes  for  o"  using  the  cross  sections 
reported  by  the  MBS  groups  (Burch,  et  al . ,  1958)  one  also  finds  that 
thissjpecies  is  apparently  in  disagreement  with  the  requirements  of 
the  "sunrise-sunset  effect."  If,  however,  the  0~2  Just  prior  to 


FIG.  IV-2  SPECTRAL  DETACHMENT  RATE  d  d  '  OF  O' AND  H' 
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detachment  is  in  a  lower-lying  state,  e.g.,  EA  =  0.46  eV,  than  that 
reported  b>  the  NBS  groups,  the  photodetachment  cross  section  could 
conceivably  be  very  small  in  the  spectral  region  above  X  =  3500  A 
(  -  3.5  eV,  and  sufficiently  large  between  2500  to  3500  A  that  radiation 
in  this  spectral  range  is  almost  entirely  responsible  for  the  photo- 
letachment  .  This  hypothesis  is  not  in  disagreement  with  Geltman’s 
theory  of  photodetachment  from  Oj. 


If  the  electron  af fimtv  of  S0~2  is  4  eV  or  perhaps  even  a  little 
larger,  photodetaching  radiation  is  restricted  to  the  near  L’V  for  this 
ion  and  it  can  meet  the  requirements  of  the  "sunrise-sunset  eifect." 
Similarly,  if  the  electron  affinities  of  Olf  and  are  -  3  eV  as  some 
measurements  seem  to  indicate,  these  ions  too  can  satisfy  the  "’sunrise- 
sunset  effect"  requirements.  Which,  if  any,  of  these  negative  ions  is 
the  dominant  one  in  the  D-rogion  can  be  ascertained  only  by  in  situ 
experiment  s . 


In  the  foregoing  we  have  considered  only  solar  radiation  in  the 
Photodetachment  process.  The  initial  burst  of  radiation  from  a  nuclear 
explosion  is  so  intense  that  the  negative  ion  concentration  will  he 
negligible  at  all  altitudes  for  several  seconds  after  detonation  due  to 
Photodetachment.  Subsequent  to  that  time  the  only  nuclear  radiations 
which  can  detach  electrons  are  and  >  rays.  For  reasonable  beta 
particle  photon  intensities  (Poppoff,  ct  id.,  1061)  ioe  to  10s  cm2  sec 
at  1  mev  energy  the  detachment  cross  sections  (essentially  the  Miller 
and^Compton  scattering  cross  sections  which  are  approximately  equal  to 
10  22  and  10-18  cm2,  respectively;  are  much  too  small  to  yield  more  than 
an  insignificant  detachment  rate.  Hence  we  conclude  that  detachment 
of  electrons  from  negative  ions  by  delayed  radiation  from  a  nuclear 

burst  is  unimportant  compared  to  solar  photodetachment  and  detachment 
by  collisions. 


Col li si onal-type  detachment  from  O2  and 
by  Whitten  and  Poppoff  (1962)  who  concluded 
polar  cap  absorption  made  by  Bailey  (1959  , 


O  was  recently  considered 
on  the  basis  of  a  study  of 
that  associative  detachment, 
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processes  ^ln/t  (lo,,  and  (lp),  is  the  dominant  mechanism.  Bailey’s 
analysis  of  PCA  events  required  a  nighttime  detachment  rate  of 
-  1CT2  sec"1.  Phelps  and  Pack  (1961)  found  that  the  rate  constant  for 
process  (1£)  was  about  4  x  10  20  cm3/ sec  at  D-region  temperatures  and 
that  molecular  nitrogen  is  ineffective  as  a  detaching  agent.  If  we 
multiply  the  rate  constant  by  the  molecular  oxygen  concentration  at  an 
altitude  of  ~  70  km  4  x  1014  cm-3),  we  obtain  a  detachment  rate  of 
only  ~  10  sec  1 .  Weber  (1962)  has  suggested  that  process  (lm)  may 
be  dominant  bv  virtue  of  the  fact  that  it  is  a  rapidly  increasing 
function  of  temperature.  However,  for  reasonable  values  of  the  tem¬ 
perature  at  this  altitude  (250-300°K)  it  is  too  small  to  agree  with  the 
empirical  value  of  the  rate  constant  suggested  by  Bailey.  Of  course,  N02 

—  *"  I 

OH  ,  or  03  may  be  the  dominant  negative  ions  in  which  case  processes 
( 1  r ) .  (Is)  or  (It)  rather  than  (lm),  ( lo )  or  (lp)  are  important. 

The  rate  constants  for  processes  (In),  (lo),  and  (lp)  have  been 
estimated  by  various  workers  to  lie  between  10_,°  and  10-15  cm1  sec-1 
vVitra,  1952;  EJalgarno,  1961).  Dalgarno’s  limiting  estimate  of 
10  1  cm3  sec  ■■1  was  presumably  made  under  the  assumption  that  almost 
every  02  -  0  and  0-0  collision  results  in  a  detaching  reaction.  This 
is,  however,  probably  too  optimistic. 

The  concentrations  of  atomic  oxygen  present  in  the  daytime  D-region 
have  been  computed  by  various  workers  (Barth,  1961;  Handbook  of  Geo¬ 
physics,  1960)  and  seem  to  lie  of  the  order  of  1011  to  10 12  cm” 3  between 
60  and  90  knu  At  n*Kht  thc  atomic  oxygen  promptly  combines  with  02  to 
form  ozone  by  means  of  a  three-body  collision  process  at  altitudes  below 
70  km.  According  to  Barth  (1961)  the  atonic  oxygen  concentration 
decreases  by  about  a  factor  of  10  at  an  altitude  of  70  km  at  night.  It 
remains  essentially  constaut  in  time  at  altitudes  of  80  km  and  above. 

If  our  effective  associative  detachment  coefficient  is  -  10-13  cm3  sec-1 
and  ,0j  10“  cm*3  at  70  km  at  night,  we  see  that  thc  nighttime  detach¬ 

ment  rate  is  -  10-2  sec'1  as  required  by  Bailey’s  analysis. 

In  addition  to  processes  ( In  , (lo  ,  and  ( lp) ,  associative  detachment 
from  ir  and  OH*  can  also  occur  (processes  ( lq  andflr  Assuming  an  atomic 
hydrogen  concentration  at  70  km  altitude  of  10‘1  cm*3  (Handbook  of  Geophysics, 
1960  ,  one  must  assign  a  rate  constant  to  the  processes  of  10'“  cm3  sec*1 
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in  order  to  obtain  agreement  with  Bailey's  analysis.  Such  values  are  not 
unreasonable.  To  our  knowledge  the  ion  0H~  has  not  heretofore  been  proposed 
as  the  dominant  negative  ion  in  the  D-region  but  the  considerations  above 
seem  to  indicate  that  it  cannot  be  excluded. 


Nighttime  detachment  from  N02  can  occur  via  process  (is)  if  the 
electron  affinity  of  this  ion  is  less  than  -  4.5  eV,  otherwise,  only 
the  particles  in  the  Maxwellian  •’tail”  of  the  velocity  distribution 
have  sufficient  energy  to  cause  electron  detachment.  Nighttime  detach¬ 
ment  from  0,  can  occur  via  process  (it).  The  values  of  the  associated 
rate  constants  are  unknown  but  kg>t  x  10“12  cms  sec'1  is  not  unreasonable 

Charge  Transfer  and  Ion-Atom  Interchange* 


In  addition  to  the  processes  of  direct  formation  and  destruction  of 
negative  ions  (la  to  lw)  there  is  the  possibility  of  charge  transfer 
and  ion-atom  interchange  (lx  to  idd).  All  of  these  processes  are  energeti¬ 
cally  possible  but  their  rates  are  completely  unknown  with  the  exception 
of  (ly).  Curran  (1962)  has  estimated  the  value  of  the  corresponding 
rate  constant  for  ( ly)  to  be  of  the  order  lO'10  to  lO'9  cm3  sec’1,  using 
the  results  of  his  mass  spectrometer  experiments.  It  is  reasonable  to 
assume  that  the  coefficients  corresponding  to  the  other  processes  (except 
lec  and  Idd  are  of  the  same  order  of  magnitude  as  the  rate  constants 
of  positive  ion  charge  transfer  reactions,  i.e.,  -  10'  > 2  to  10"°  cm2  SCc-. 
Processes  (lcc  and  idd)  involve  the  breaking  and,  or  forming  of  chemical 
bonds  and  may  thus  have  considerable  activation  energies.  If  so,  they 
are  expected  to  be  slower  than  the  charge  transfer  reactions.  Actually, 
process  lcc)  is  energetically  forbidden  unless  NO^  has  an  electron 
affinity  greater  than  „  5.7  eV .  However,  it  has  been  previously  shown 
that  nighttime  detachment  will  not  occur  unless  the  electron  affinity 
is  less  than  ~  4.5  eV.  We  therefore  rule  out  reaction  (lcc) 


.Similar  processes  (lee  and  iff)  can  produce  0'3  ; 
Since  molecular  oxygen  is  much  more  effective  than 
in  three-body  attachment  (la),  it  is  probable  that 
consideration  of  third  bodies  to  02 . 


both  are  exothermic, 
molecular  nitrogen 
we  can  limit  our 
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Negative  Ion  Kinetics 


5  . 


In  the  foregoing  we  concluded,  for  various  reasons,  that  many  of 
the  negative  ion  processes  can  probably  be  neglected.  In  addition,  it 
can  be  shown  that  the  ion-ion  recombination  processes  above  80  km  are 
probably  much  slower  than  the  others  listed.  Those  remaining  are  (la), 
(le).  (IS),  (lh),  (lj),  (lk),  (In)  to  (lp),  (lr)  to  (It),  (lx),  (ly), 
(laa),  (lbb),  and  (ldd  to  (iff);  the  kinetic  equations  for  these 
processes  can  be  written: 

d 

dt*  J  =  K  °2  2  0  -  f(°2.  02  ~  >  02]  0]  -  k  [ojKo]  -  k  [ 02 ] [ OH ] 

n  >  °  x  aa 

“  kddiO;j[N][N2] 


-  kff [02][0][02]  (12) 

~  =  a:o][cj  +  kx  ojlo]  -  o(o')l0"  -  7p  o"][0]  -  k  [0“)[N0,  ] 


dtN02. 

dt 


d,0H~ 

dt 

d  i  Oa  j 
dt 


-  koe  °'][o2;2  (i3) 

ky‘°  N°2  f  kdd’°*1,N  N'2  -  0(N0^)(N07j  -  7B(N0aJ(N]  (14) 

kaa  °2  0,1  +  kbb  0  0,1  "  L(0,O  0,i~]  -  >r  Oifj  H  (15) 

kee  0  °2jl  +  kff  °2  0  °2J  “  n(02)  0 3  J  -  >t[o2  0  (16) 


where  K  is  the  rate  constant  for  three-body  electron  attachment  to  02 , 
p  is  the  coefficient  of  radiative  attachment  to  0,  c  is  the  photo¬ 
detachment  rate  for  the  species  indicated,  e.g.,  o(02),  y  are  the 
coll l si onal  type  detachment  rates,  and  k  the  rate  constants  for  the 
charge  transfer  processes;  e.g.,  ky  is  the  rate  constant  for  process 

(lv),  etc.  Square  brackets  represent  the  concentrations  of  the  various 
species . 
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The  following  values  of  some  of  the  rate  parameters  are  Assumed; 

K  =  1.5  x  10  30  cm3/sec  (Chanin,  et  1959) 

p(02)  =  o.l  sec-1  (this  paper) 
p(0  )  =  1.4  sec-1  (Branscomb,  et  al,,  1958) 
p(N02)  =  10  6  sec-1  (assumed) 

P  =  1.34  x  10" 15  cm3  sec-1  (Branscomb,  et  aK ,  1958) 


kv  ~  3  x  10'10  cm3  sec-1  (Curran,  1962) 


**dd  kee  *  kff  =  10  28  cmS  sec-1  (assumed) 

In  addition  we  assume  daytime  concentrations  of  N2 ,  02 ,  0,  N02 
OH,  H  and  electrons  at  75  km  altitude  as  follows; 

N2]  =  8  x  10* *  cm-3  (Minzner,  1959) 

02 j  =  2  x  101 *  cm-3  (Minzner,  1959) 

O  =  2  *  10“  cm-3  (Barth,  1961) 

N02  1=  10J  cm-3  (Barth,  1961) 

OH  1°  cm  3  (Handbook  of  Geophysics,  I960) 

H  10  cm  (Handbook  of  Geophysics,  1960) 

e.  =  500  cm-3  (day)  (Galejs,  1962) 

100  cm"3  (night)  (Galejs,  1962) 


—  . negative  ions  tor 

P  ble  values  of  ?t ,  t  (OH  ),  p(03),  and  tN]  are  presented  in  Table  jv-1 
Although  the  set  of  values  used  is  by  no  means  exhaustive,  it  does  give 
some  idea  of  the  range  of  negative  ion  concentrations  which  may  be 
expected  under  daytime  and  nighttime  conditions.  It  is  apparent  from 

the  right-hand  columns  that  OH',  NO^ ,  and  0'3  may  be  of  great  importance 
in  the  D-region. 
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The  var  ous  rate  constants  suggested  for  Off,  0^ ,  and  NoJ  are 
entirely  reasonable.  For  example,  the  photodetachment  rate  for  OH- 
could  well  be  10-  sec-  if  the  electron  affinity  is  as  large  as  3  eV 
It  is  quite  possible  that  all  four  species  o'.  Off,  0;,  and  NO^  are 
important  in  the  D-region,  the  relative  importance  of  each  varying 
with  altitude  and  condition  of  daytime  or  nighttime.  We  re-emphasize 
that  the  foregoing  are  only  plausibility  arguments. 


6.  Conclusions 


The 


following  conclusions  are  drawn  from  the 


preceding  arguments. 


a . 


One  cannot  identify  the  species  o  and  ff 
negative  ions  in  the  lower  D-region,  since 
with  the  "sunrise-sunset  effect."  o'  may 
the  upper  D-region. 


with  the  dominant 
they  both  conflict 
be  important  in 


b. 


02,  OH 


03,  or  NOj  may  be  the  dominant  species  of  negative 
ion  without  contradicting  existing  knowledge  of  the  D-region. 
In  fact  they  may  all  be  important. 


c .  Radiative  attachment  to  02  and  three-body  attachment  to  O 
are  probably  not  important  processes.  Further  investigation 

is  required  before  this  statement  can  be  definitely  established, 
however. 

d.  Associative  detachment  is  probably  the  dominant  nighttime 
mode  of  electron  detachment  from  the  species  0~  and  OH 
Detachment  from  N02  and  03  can  occur  via  dissociative 
processes  ^ls  and  It). 

».  Previously  proposed  values  of  the  photodetachment  rate  for 
02  are  not  reliable. 

f.  Diffusion  to  dust  i s  an  unimportant  mode  of  negative  ion 
removal  (see  Section  VIII)  . 

At  the  risk  of  seeming  trite  we  reiterate  the  need  for  further 
measurements,  both  laboratory  and  in  situ,  if  a  reliable  model  of 
negative  ion  species  and  distribution  is  to  be  constructed. 
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B‘  --arge  Transfcr  and  Ion-Atom  Interchange  (Positive  Ions) 

In  this  section  we  shall  investigate  the  importance  of  charge  transfer 
and  ion-atom  interchange  processes  with  respect  to  the  ionic  and  electronic 
structures  of  the  D-  and  E-regions.  Specifically  we  shall  show  that  these 
processes  control  the  concentrations  of  the  various  ion  species  by 
’’shuffling"  them  more  rapidly  than  they  can  be  created.  The  nature  of 
the  ionic  species  in  turn  has  an  important  bearing  on  the  rate  of  electron 
recombination  and  hence  the  electron  density  profile.  There  are  many 
such  processes  which  are  possible  in  the  altitude  interval  80  to  150  km 
.Bortner,  1961a;  Bortner  and  Baulknight,  1961b;  Nawrocki .  1961a)  but 
-st  can  be  eliminated  for  various  reasons,  e.g.,  a  participating  neutral 
species  is  present  only  in  trace  concentrations,  reactions  are  highly 
endothermic,  activation  energies  are  too  large  for  the  process  to  be 
important  at  ordinary  temperatures,  etc.  The  following  reactions  are 
considered  potentially  important  in  the  shuffling  process. 


+ 


+  o2 

0 

+ 

0 

+ 

(1) 

+  0 

-  n2  +  0+ 

(2) 

4  02 

•  N2  +  02 

(3) 

+  NO 

-  N2  4  N0+ 

(•1) 

+  NO 

-  02  4  N0+ 

(5) 

+  NO 

-  0  4  N0+ 

(6) 

4  n2 

•  NO*  t  K 

(7) 

+  NO 

-  02  4  N 

(«) 

4  N 

-  NO*  4  0 

0) 

+  02 

-  NO  4  N0+ 

(10) 

+  n2 

•  NO  4  N0+ 

(II) 

+  0 

•  N  4  N0+ 

(12) 
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However,  two  of  the  processes  [(10)  and  (ll)'J  involve  the  breaking 
of  two  double  bonds  and  hence  are  expected  to  have  large  activation 
energies.  Bortner  (1961a)  has  assigned  rate  constants  of  7  x  10~,s  and 
3  x  10-17  era3  sec,  respectively,  but  these  are  admittedly  crude  estimates 
only  and  are  probably  upper  bounds  on  the  actual  values.  In  any  event 
these  estimates  arc  sufficiently  small  that  processes  no)  and  (11)  can 
be  neglected  in  any  consideration  of  the  ionosphere  above  90  km.  Fite 
and  coworkers  (Fite,  Rutherford,  Snow  and  van  Lint,  1962)  have  found 
that  the  former  is  much  slower  than  the  charge  transfer  process  (5). 

i the  less  reaction  Uj  must  be  considered  in  ion  shuffling  kinetics 
at  altitudes  below  90  km.  D.  R.  Bates  (private  communication)  has 
concluded  that  process  (12)  is  suppressed  bv  virtue  of  the  fact  that 
an  electronic  transition  is  involved.  have  therefore  assumed  it  to 
be  of  negligible  importance  in  our  model. 


The  largest  possible  value  of  the  rate  constants  for  any  of  the 
processes  (1)  to  (9)  occurs  if  a  reaction  takes  place  during  every 
collision;  this  value  is  of  the  order  or  10“ 9  cm3/sec  if  we  equate  the 
react  ion  cross  section  to  the  classical  collision  cross  section  tyring, 
Hi rsh folder,  and  Taylor,  1936).  For  reasonable  estimates  of  the  concen¬ 
tration  of  nitric  oxide  in  the  upper  atmosphere  (see,  for  example, 
Handbook  of  Geophysics,  I960)  the  reaction  rates  of  processes  j.|),  (5), 

l'iJ,  iU,<l  :,r‘*  to°  snMil  these  processes  to  be  important  in 

determining  the  equilibrium  concentrations  of  the  principal  i„„  species. 
Process  v9)  is  believed  unimportant  for  a  similar  reason  (small  Nj). 

The  remaining  processes  are  (1  ,  (2).  (3),  (7)  and  below  90  km,  (n), 

:U"'  tlu‘  01,11  *  librium  equations  for  the  various  ion  species  can  be 
approximated  by 


“a 

a 

U|) 

N\.  *  J  « 

j  i  k2  N2  *  j  ;  o 

+  k  3  N2  0  j 

(13) 

b 

•  .  -*  f 

4- 

“b 

‘l> 

().,  , • 

.  -  k3  N.,  j  o2 

rq 

o 

1 

-  ki,  o2'j  \2J 

k7 

0+  j[N2J 

'  k,  0+j  02 

-  k  2  f  N  2  '  j  '  0  j 

US) 

0  = 

d 

‘d 

NO  c 

0  j  [  N2 

kjj  02  N2 

(16) 
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where  qfl  c  represent  the  production  rates  of  N2 ,  0  2 ,  and  0+,  respec¬ 
tively,  the  k 's  are  the  rate  constants  corresponding  to  certain  of 
the  processes  (as  indicated  by  the  subscripts)  (l)  to  (9),  the  square 
brackets  represent  the  concentrations  of  the  indicated  ions  and  neutral 
species,  and  ’  ’  represent  the  dissociative  recombination  coefficients 
of  N2 ,  02 ,  and  NO  ,  respectively.  We  shall  show  in  section  iv-D  that  ion- 
ion  recombination  can  probably  be  neglected  in  the  daytime  at  altitudes 
above  80  km.  Obviously  we  must  require  that 

[ej  =  N2  J  +  [02]  +  [0+j  +  [NO+]  (17) 

and 

q  (electron)  -  %  +  %  +  Q  =  q  (ion)  (18) 

if  electric  charge  is  to  be  conserved. 


In  addition  to  the  effects  of  recombination,  two  other  electron- 
ion  removal  processes,  diffusion  and  ionospheric  currents,  should  in 
principle  be  included  in  the  equilibrium  equations.  The  effects  of 
vertical  diffusion  were  computed  by  finding  the  approximate  value  of  the 
term  r^(N+W)  where  is  the  positive  ion  density,  h  refers  to  altitude 
and  W  is  the  diffusion  velocity  given  by 


W  =  -D(i  ™  +  I) 

vn  >h  ir 


(19) 


In  equation  (19  D  is  the  diffusion  coefficient,  N  is  the  neutral  species 
concentration  and  His  scale  height.  According  to  Nawrocki  and  Papa 
(19(>lb;  the  diffusion  coefficient  of  NO  in  air  (at  temperatures  >  300°K) 
is  given  approximately  by 

-J 


D  = 


9  x  1016  T2 
N 


(20) 


where  T  is  the  temperature  in  °K .  Using  reasonable  values  of  T  and  N 
.Minzner,  e£  al. ,  1959  ,  it  was  found  that  the  diffusion  term  was  only 
about  four  percent  of  the  magnitude  of  the  recombination  term  at  an 
altitude  of  150  km.  Hence  the  former  1=  entirely  neglected  in  this 
treatment  since  it  rapidly  decreases  as  one  descends  to  lower  altitudes. 
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The  effects  of  the  ionospheric  current  system  on  electron-ion 
removal  have  only  recently  been  discovered  (L.  H.  Brace,  private  communi¬ 
cation).  At  the  present  time  we  do  not  have  sufficient  information  for 
a  quantitative  study  of  its  effects  on  the  ion  concentrations  present 
in  the  E-region.  Since  the  rate  constants  reported  here  are  order  of 
magnitude  only,  it  is  not  likely  that  they  would  be  greatly  changed  by 
inclusion  of  the  current  term.  This  point  will  be  investigated  in  a 
future  report  when  more  data  are  available. 

If  the  values  of  the  various  rate  constants  as  well  as  the  q's 
are  known,  the  ion  concentrations  can  be  readily  computed.  Conversely 
the  rate  constants  k  can  be  found  if  the  ion  concentrations,  ion  production 
rates,  and  dissociative  recombination  coefficients  are  known .  Actually, 
sufficient  knowledge  of  the  upper  atmosphere  for  the  determination 
of  precise  values  of  any  of  these  quantities  does  not  exist.  Never¬ 
theless  one  can  obtain  order  of  magnitude  estimates  of  the  rate  constants 
from  mass  spectrometric  measurements  of  the  relative  concentrations  of 
various  ionic  species  between  100  and  200  km  altitude  together  with 
estimates  of  the  dissociative  recombination  coefficients  aa,b,d  based  on 
results  obtained  in  the  next  section  and  the  production  raLs  of  the 

ionic  species,  The  mass  spectrometric  measurements  of 

the  relative  concentrations  of  the  various  species  of  positive  ions  at 
various  altitudes  (Johnson,  Meadows,  and  Holmes,  1958;  Meadows  and  Townsend 
I960;  Taylor  and  Brinton,  1961;  Istomin  and  Pokhunkov,  1962)  are  listed 
m  Table  IV-2.  The  values  of  the  dissociative  recombination  coefficients 
presented  in  Table  IV-3  are  based  on  considerations  discussed  in  the  next 
section.  The  concentrations  of  N2 ,  02 ,  and  0  at  various  altitudes  accord¬ 
ing  to  a  model  developed  at  the  National  Bureau  of  Standards  (Norton,  van 
Zandt  and  Denison,  1962)  arc  presented  in  Table  IV- 1.  Production  rates 
<>f  N2,  02,  and  0+  Norton,  van  Zandt,  and  Denison,  1962)  based  on  recently 
obtained  solar  ultraviolet  spectra  (Watanabe  and  Hinteregger,  1962)  are 
listed  in  Table  IV-5 
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Table  IV- 2 


RELATIVE  CONCENTRATIONS  OF  POSITIVE  ION  SPECIES 


Rocket 


Fraction  of  Total  Positive  Ions  which  ai 


Flight 

(km) 

.+ 

n2 

ot 

N0+ 

0+ 

NN  3.17  F1 

150 

<0.03 

0.-1 

0.5 

0. 1 

2321  CST 

120 

<0.03 

0.2 

0.8 

<.0.03 

20  Nov.  56 

100 

-  0.03 

-  0 

1  .0 

-<  0 

NN  3.18  F1 

150 

0.03 

0.1 

0.8 

0.1 

2002  CST 

120 

0.03 

0.1 

0.9 

0.03 

21  Feb.  58 

100 

0.03 

0.1 

0.9 

•<  0 

NN  3.19  F* 

150 

0.03 

0.3 

0.6 

0 . 1 

1207  CST 

120 

<0.03 

0.3 

0.7 

0.03 

23  Mar.  58 

100 

0.03 

0.1 

0.6 

0 

NASA  -1.09 2 

150 

0.01 

0.3 

0.6 

0.1 

10-17  EST 

120 

0.01 

0.1 

0.9 

0.05 

20  Apr.  60 

100 

0.01 

0.  1 

0.9 

0 

NASA  -4.1-12 

150 

0.01 

0.-1 

0.1 

0 . 2 

11-11  EST 

120 

0.01 

0.2 

0.8 

0.05 

15  Nov.  60 

100 

0.01 

-  0 

1.0 

0 

I  S  SR  3 
.  .o 

150 

— 

0.1 

_  _ 

h  =  —  b  o 

120 

-- 

_  _ 

9  Sept.  57 

100 

— 

— 

— 

_  _ 

ISSR  3 

150 

— 

0.2 

0.6 

0.2 

h  =  36  m 

120 

-- 

0.3 

0.7 

o 

2  Aug.  58 

100 

— 

0.2 

0.8 

0 

USSR  3 

.  O 

150 

— 

0.1 

0.6 

0.05 

h  =  0  m 

120 

— 

0.2 

0.8 

^  0 

13  Aug.  58 

100 

— 

-- 

— 

ISSR3 

150 

-- 

0.3 

0.7  ' 

0.05 

h  =  0  m 

120 

_ 

0.1 

0.9 

o 

22  July  59 

100 

-- 

— 

ISSR3 

150 

— 

0.-1 

0.5 

0. 1 

h  =  15  m 

120 

— 

0.3 

0.7 

0 

15  June  60 

100 

0.1 

0.9 

-  0 

‘Johnson,  Meadows,  and  Holmes,  1958;  Meadows  and  Townsend,  1960 
2Tay lor  and  Brinton,  1961 


Table  IV-3 


DISSOCIATIVE  RECOMBINATION  COEFFICIENTS 

(cm3  sec-1) 


Table  IV- 5 


ION  PRODUCTION  RATES3 
(cm'3  sec-1) 


I  Altitude 
(km) 

02 

0 

200 

*150 

-  0 

800 

150 

2000 

160 

1700 

120 

1200 

600 

800 

100 

■100 

1100 

160 

h 

80 

50 

50 

-  0 

a 

Norton,  van  2andt,  and  Denison,  1962 


Estimated  ionization  rates  under 
disturbed  solar  conditions. 


In  order  to  arrive  at  estimates  of  the  rate  constants  k,,  k,  ,  k3, 
and  k7  it  was  first  necessary  to  compute  the  electron  concentration  at 
the  altitudes  considered.  This  was  done  by  setting  the  sum  of  the  right- 
hand  members  of  equations  13;  to  (16)  equal  to  the  total  ion  production 
rate,  expressing  the  various  ion  concentrations  in  terms  of  the  electron 
concentration  ffrom  Table  IV-2)  and  solving  for  e,.  The  rocket  data 
employed  in  the  computations  were  taken  from  Nights  NN  3.19  (solar 
zenith  angle:  i  =  32°),  NASA  1.09  (X  =  28°),  and  the  Soviet  flight  of 
2  August  1958  (x  =  5i°).  Although  the  solar  zenith  angle  in  these  flights 
was  not  equal  to  the  ,  =  3-1°  for  which  the  ion  production  rates  were  com¬ 
puted,  this  discrepancy  introduces  very  little  error,  particularly  in 
the  case  of  the  American  observations.  The  solar  zenith  angles  of  the 
other  observations  were  so  large  that  one  could  not  justify  their  use. 

Values  of  the  rate  constants  at  an  altitude  of  150  km  were  computed 
bv  substituting  the  values  of  the  estimated  ion  and  electron  concentra¬ 
tions  and  the  ion  production  rates  from  Table  1V-5  together  with  the 
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appropriate  dissociative  recombination  coefficients  (Table  IV-3)  and 
neutral  species  concentrations  (Table  IV-4)  into  equations  (13)  to  (16). 
The  values  which  best  fit  the  data  were  found  to  be 


k, 

-  2 

X 

1 

o 

H 

cm3  see" 

k2 

2 

X 

10“'  1 

cm3  sec¬ 

k  3 

2 

X 

10“i° 

ern3  sec“ 

kr 

6 

X 

n 

1 

O 

cm3  sec"1 

Of  those,  k , .  k3,  and  k7  have  recently  been  measured  in  the  laboratory 
by  various  workers.  Th<  results  are 


k,  nickinson  and  Sayers  (i960) 
Sayers  (1962) 

I.angstroth  and  Hasted  (1962) 

Fite,  Rutherford,  Snow,  and 
van  Lint  (1962) 

k 3  Fite,  Rutherford,  Snow,  and 
van  l.i  nt  (1962) 

k7  Tal rose ,  Markin,  and  Larin 
(1962 

Sayers  V1962; 


2.5  x  10“ 1 1  cm3  sec 

1.6  x  10“' 1 
1.6  x  10“ 12 

(1  to  15)  x  10“* 1 

2  x  10“ 10 

6.7  x  10“»2 
2.5  x  10-' 1 


I.angstroth  and  Hasted  (1962) 


•1.7  x  10“12 


It  is  evident 
agreement  with 


that  our  estimates  of  k,  and  k3  are  in  quite  good 
the  laboratory  measurements.  In  fact  the  agreement 


i  s 


,  r-  -  ....  v-nuillf. 


. h"*  veil  UL'fi 

only.  On  the  other  hand,  our  estimate  of  k  ,  , 

,  ou>  esumate  of  k7  is  in  very  poor  agreement 

w  1,1  tlle  laboratory  measurements.  Since  the  nstimt  h 

oiine  mi  estimated  uncertainty  in 


our  computations 
conclusions  are  ( 
using  too  small  a 


is  only  about  half  an  order  of  magnitude,  the  obvious 
tf  the  laboratory  measurements  are  correct)  that  we  an 
value  for  the  dissociative  recombination  coefficient 


74 


of  NO  or  that  another  charge  transfer  process  for  conversion  of  N0+  into 
a  different  ion  is  effective.  However,  it  is  believed  highly  improbable 
that  Q^(NO+)  is  in  error  by  more  than  half  an  order  of  magnitude.  Nor 
can  we  find  an  exothermic  charge  transfer  process  which  will  remove  the 
NOj  all  such  exothermic  processes  involving  NO+  art  ones  which  produce 
NO  .  Hence,  the  discrepancy  remains  unresolved. 


The  concentrations  of  .<,  02,  0+,  and  NO+  at  various  alti 
been  computed  using  our  estimated  rate  constants  and  the  means 
laboratory-determined  values  listed  above.  They  are  shown  in 


tudes  have 
of  the 
Figures  IV-3 


riv.  iv-j 


ION  CONCENTRATIONS  ABOVE  80  km.  THESE  CONCENTRATIONS  WERE 

COMPUTED  USING  THE  VALUES  OF  THE  RATE  CONSTANTS  k  k  k 

AND  k7  OBTAINED  FROM  THE  MASS  SPECTROMET RIC  OBSERVATIONS 
THE  BROKEN  LINE  REPRESENTS  THE  OBSERVED  CONCENTRATION  OF  NOf 


and  T  V— 1  .respectively .  The  broken  line  in  Figure  IV-3,  which  v.as  ob¬ 
tained  from  the  observed  relative  concentration  of  NO+ ,  illustrates  the 
discrepancy  between  the  computed  ion  concentrations  using  our  set  of  rate 
constants  and  the  observed  relative  ion  concentrations  at  altitudes 
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H  3-  4  M 


FIG.  IV-4  ION  CONCENTRATIONS  ABOVE  80  km.  THESE  CONCENTRATIONS  WERE 

COMPUTED  USING  THE  MEANS  OF  RECENT  LABORATORY  MEASUREMENTS 
OF  k,,  k3,  AND  k7.  OUR  VALUE  OF  k  WAS  USED  IN  THE  COMPUTATIONS 


UMo,  150  km.  Comparison  of  the  NO *  curve  ,n  Figure  IV-I  with  tin  observed 
relativ,  NO  *  concentrations  reveals  the  d  ,  sagreemen  l  between  observation 

and  our  model  when  t  he  rate  constant  k7  and  Savers  (1962)  is  used  in  the 
computations . 


It  is  evident  from  Figures  IV-3  and  IV-I  that  tin  dominant  ion  in 
the  D-rcgion  under  disturbed  conditions  must  be  o’.  NO*  is  quite  unim¬ 
portant  unless  it  is  produced  directly,  e.g.,  ionization  of  NO  by  i.yman  a. 
This  model  cannot  be  extended  to  the  lower  D-rcgion  (-'  75  km)  because  of 
the  possible  importance  of  such  exotic  ions  as  N, ,  N*  ,  and  0.;.  Kasner, 

—  — •  11961  havc  obsL'rvt’d  such  ions  at  gas  pressures  of  about  0.1  mm  of 
HB-  0bla,ni,,B  3  va,ue  ~  ^  x  10-e  cm3  sec-  at  this  pressure. 

However,  the  pressure  at  which  such  ions  gain  significance  as  well  as  the 
c  fects  at  higher  pressures  (e.g.,  formation  of  etc.)  are  not  known. 
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C.  Dissociative  Recombination 


Numerous  laboratory  measurements  (see  TableIV-6)  of  the  dissocia¬ 
tive  recombination  coefficients  of  N2  and  02  have  been  made  during  the 
last  decade  with  somewhat  conflicting  results.  In  addition,  three 
measurements  of  the  recombination  coefficient  of  N0+  have  been  carried 
out  recently;  one  of  these  was  a  shock  tube  measurement  (at  a  tempera¬ 
ture  of  5000°K ) ,  and  the  others  were  low  temperature,  i.e.,  room 
temperature,  plasma  experiments.  Since  most  of  these  laboratory  values 
were  discussed  earlier  by  Bortner  and  Baulknight  (1961),  the  purpose 
of  this  report  is  to  include  some  new  results  (Kasncr,  et  al . ,  1961; 

Cunton  and  Inn,  1961;  Doering  and  Mahan,  1962)  together  with  the  results 
of  various  ionosphric  measurements  (McElhinny,  1959;  Bowhill,  1961;  Whitten 
and  Poppoff,  1961)  and  to  develop  from  them  a  model  of  the  daytime  ef fee- 

recombination  coefficients  in  the  altitude  interval  80  to  150  km. 

The  recent  measurement  of  cyN*)  and  cyo2)  by  Biondfs  group 
involved  the  use  of  a  mass  spectrometer.  They  found  that  at  low  gas 
pressures,  i.e.,  0.01  mm  Hg  of  h’2  and  02  ,  the  dominant  positive  ions 

were  N2  and  02 ,  but  at  higher  pressures,  i.e.,  >  0.1  mm  Hg  of  N2  and 
02,  the  dominant  positive  ions  were  it*,  .<,  and  0+, .  This  was  in  good 
agreement  with  the  earlier  results  of  mobility  experiments  conducted  by 
Saporaschenko  (1958,  and  Kovar  and  coworkers  (1957).  The  latter  investi¬ 
gators  also  found  a  strong  temperature  dependence  for  the  formation  of 
N«:  its  formation  is  greatly  inhibited  by  high  temperatures.  This  is 

readily  understandable  in  view  of  the  low  binding  energy  of  ,  about 
0.1-1  eV. 


Measurements  of  the  decrease  in  the  E-region  critical  frequency 
during  colar  eclipses  have  been  used  by  McElhinny  (1959)  to  obtain 
effective  daytime  recombination  coefficients  at  E-region  altitudes. 

Bow  v19^1  has,  by  use  of  a  two-ion  model  and  McElhinny's  effective 
coefficients,  derived  component  recombination  coefficients.  The  values 
were  a ‘  =  10“?  cm3  sec"1  and  a3  =  5  x  10"9,  ^  =  1  x  10"*  cm3  sec’*, 
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Table  IV-6 


LABORATORY  MEASUREMENTS  OF  DISSOCIATIVE 
RECOMBINATION  COEFFICIENTS  Q 


Investigators 

VNa) 

(cm3  sec-1 ) 

|  ^q(°2  ) 

(cm3  sec-1) 

qd(no+) 

(cm3  sec-1) 

Kasner,  Rogers,  and  Blondi ,  1961 

a(5 . 9+ 1 )  x  lO"7 

a(3.8±l)  x  10-7 

b 

2  x  10-6 

Holt,  1959 

3  x  10~7 

Falre  and  Champion,  1959 

4.0  x  10~7 

Blondi  and  Brown,  1949 

1,4  x  10“e 

Falre,  Fundingsland,  Aden,  and 
Champion,  1958 

1.4  x  10-s 

Bialecki  and  Dougal ,  1958 

C6.7  x  10“6 

d8.7  x  lO"7 

Sayers,  1956 

Lin,  I960 

Gunton  and  Inn,  1961 

Doering  and  Mahan,  1962 

1 

*1.1  x  10~7 

f4  x  10~8 

J( 

Bio-* 

h1.3  x  10~* 

3^?)  x  lO"7 
—  _ 

afor  pressures  <  O.Ol  mm  Hr;  temp  -  300°K 
’for  pressures  >  0.1  mm  Hr;  temp  -  300°K 


c  O 
at  92  K 

d  q 

at  300  K 
e  o 

at  3200  K 

f  o 

at  2500  K 

Kat  500°K 

pressures  >0.1  mm  Hg;  temp  -  300°K 
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presumably  representing  recombination  of  02  and  N0+,  respectively; 
the  two  values  of  Ct3  correspond  to  different  assumptions  as  to  the 
altitude  of  emission  of  ionizing  radiation  from  the  solar  corona.  In 
his (1961)  paper  Bowhill  also  computed  values  of  and  a3  using  measure¬ 
ments  of  the  diurnal  variation  of  the  E-region  critical  frequency  together 
with  his  two-ion  model.  The  component  recombination  coefficients  were 
found  to  be  =  10-7  cm3  sec"1  and  a3  =  6  x  10-9  cm3  sec'1.  This  value 
of  a3  tends  to  support  the  first  one  listed  above  for  eclipse  observations. 

The  computation  of  an  effective  dissociative  recombination  coefficient 
for  the  D-region  during  sudden  Ionospheric  disturbances  was  presented  in 
III-A;  =  2t3  x  10  7  cm3  sec-1.  It  was  shown  in  the  preceding 

section  that  we  do  not  know  whether  the  dominant  ion  in  the  disturbed  D- 
region  is  02  or  NO  .  Hence  we  do  not  know  what  ratio  of  [Oj],  [N0+’  is 
represented  by  our  value  of  the  recombination  coefficient. 


Dissociative  recombination  coefficients  are  temperature-dependent, 
decreasing  with  increasing  temperatures.  Several  years  ago  Bates  (1950) 
derived  the  well-known  form  for  general  temperature  dependence,  a  cc  t'^. 
Recently  a  more  general  form  was  derived  under  the  assumption  of  thermal 
equilibrium  (Bates  and  Dalgarno,  1962)  which  indicates  that  the  tempera¬ 
ture  dependence  of  ^  la  mjch  more  compUcatc()  As  thpy  pol[U  out  >  ^ 

transition  may  occur  with  or  without  emission  of  radiation,  the  latter 

case  being  the  more  common.  Their  form  of  the  dissociative  recombination 
coefficient  is 


*D  -  Cl 


u> 

3/  Z  AD 

00.  4 

AB 


Ao  (T) 


(1) 


where  u)^  and  u>ab+  are  the  statistical  weights  of  the 
states  of  the  neutral  molecule  and  the  molecular  ion, 


relevant  electronic 
and  C,  is  a  constant. 


C, 


h3 

2(2-mk  3/3 


(2) 


in  which  h,  k,  and  m  are  Planck's  and 
of  the  electron,  respectively.  A0(Tl 


Boltzmann's  constants  and  the  mass 
is  a  temperature-dependent  function 
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related  to  the  oscillator  strength  of  the  transition,  f0(e),  and  to  the 
lifetimes  of  the  radiation  emission  and  radiationless  modes  (r  and  r  ) 
by  the  equation  P  a 


where 


A0(T) 


_ go(T) 

Ta  +  Tp  go(T) 


(3) 


Ko(T)  =  /  f o  (e  )  e  kT  de 
o 


(4) 


The  subscript  zero  affixed  to  the  oscillator  strength  indicates  that  we 
are  considering  the  system  to  be  in  the  ground  vibrational  state.  If 
the  temperature  is  high,  excited  vibrational  states  v  must  be  considered; 
the  function  g0  must  then  be  replaced  by  an  average 


K  (T) 


£  ev(T) 


(5) 


For  E-region  temperatures,  excited  vibrational  states  cannot  always  be 
neglected. 


nccoramg  to  Hates  and  Dalgarno,  the  case  in  which  r  g0(T)  «  t 


is  the  more  common.  Then  A0  takes  the  form 


Ao(T)  =  —  j  f0(t)  e  t//kT  de 
a  o 


(6) 


The  temperature  dependence  of  0<D  is  thus  dependent  in  turn  upon  the  func¬ 
tional  form  of  f0(e).  For  example,  if  the  latter  is-e~*  ,  then  1  cc  T" 1 
fn  general,  however,  one  would  expect  the  energy  dependence  of  f0(e)  to 
be  more  complex  than  a  power  law  in  which  case  the  dissociative  recombina¬ 
tion  coefficient  would  not  be  related  to  the  temperature  by  any  simple 
expression,  particularly  at  elevated  temperatures.  Nevertheless,  values 
of  the  exponent  n  in  the  range  -3,  2  <  n  <  -1,2  seem  to  fit  the  existing 
experimental  data  fairly  well.  This  is  shown  in  Figures  IV-5,  6  and  7 
which  indicate  that  n  .  -3,4  for  0^)  and  n  -  -1  for  ^(oj)  are  the 
best  fits  to  the  most  reliable  data.  Actually  the  temperature  exponent 
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,.3»»/jlUD - (  JO  )°C 


T  f*K)-' 


FIG.  IV-5  DISSOCIATIVE  RECOMBINATION  COEFFICIENT  OF  N*  AT  VARIOUS 
TEMPERATURES  T  2 


T*  — (-K)  ' 


FIG.  I V-6  DISSOCIATIVE  RECOMBINATION  COEFFICIENT  OF  0*  AT  VARIOUS 
TEMPERATURES T  2 
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corresponding  to  Nj  may  be  as  small  as  -1  or  as  large  as  -1  2.  In 
view  of  Kasner’s  results  it  would  appear  that  the  poor  agreement  of  the 
data  of  Bialecki  and  Dougal  (1958)  is  due  to  the  formation  of  relatively 
laige  concentrations  of  N3  and  N4  in  their  experiments.  Figure  IV-G 
also  presents  the  previously  mentioned  values  of  cy  presumably  corre¬ 
sponding  to  recombination  with  ot .  obtained  from  eclipse  and  diurnal 
variation  observations  (McEl hinny,  1959;  Bowhill,  1961).  The  temperature 
to  Which  this  value  corresponds  can  only  be  estimated  but  the  agreement 
with  laboratory  measurements  is  quite  good 

Derivation  of  the  temperature  dependence  of  a  (N0+)  is  much  more 
tenuous  and  is  based  on  the  value  obtained  by  I. in  (1960)  and  Doering 
and  Mahan  (1962)  as  well  as  on  solar  eclipse  data  (McElhinnv,  1959; 
llowhill,  1961).  An  approximate  value  of  n  Tor  u()(N0+)  is  n  -3,2. 

Again  the  temperatures  to  which  the  values  obtained  from  ionospheric 
observations  correspond  must  be  estimated  (~  1()00°K  Tor  the  K-region) 

Ciunton  and  Inn  (1961)  have  obtained  a  value  for  cyNO*)  or  1.3  x  10'6  cm3  sec 
at  room  temperature  and  at  nitric  oxide  pressures  of  several  tenths  or 
one  mm  of  llg  This  value  seems  too  large  to  correspond  to  the  dissociative 
recombination  or  No"*  and  may  correspond  to  recombination  or  complex  ions 
such  as  NO  2  Somewhat  similar  results  were  obtained  by  Dooring  and 
Mahan  who  estimated  the  experimental  upper  limit  as  2  x  ]()~B  cm3  sec 
and  the  lower  limit  as  10’7  cm3  sec.  The  mine  ox,de  partial  pressure 
was  of  the  order  of  0.1  mm  of  llg. 

The  conclusions  concerning  dissociative  recombination  can  be 
summarized  as  follows; 


1 


On  the  basis  of  existing  data,  probable  values  of  the  impor¬ 
tant  dissociative  recombination  coefficients  in  the  D-  and  E- 
regions  arc  y  N* )  =  5  x  10"5  1" 3  4  cm3  sec-1;  ^  o\  )  = 


9  x  10“5  T-'  cm3  sec"1;  cyNo\  =  2  6  x  10 


3  2  cm3  sec- 
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2.  It  is  evident  from  the  results  of  the  preceding  section  that 
the  most  important  coefficients  in  the  upper  D- region  and 
in  the  E-region  are  cyctj)  and  cyNO4);  is  at  most  of 

marginal  significance. 

3  The  identities  of  the  dominant  positive  ions  in  the  lower  D- 
tegion  are  not  known.  This  arises  from  ignorance  of  the  role 
played  by  the  ions  ,  Nj,  and  0,  at  these  altitudes.  (See 
Section  IV-B,  ion-atom  interchange  processes.) 

A  model  of  the  effective  daytime  ionospheric  electron  dissociative 
recombination  coefficient  in  the  altitude  range  80  to  150  km  is  shown 
in  Figure  IV-8  The  profiles  A,  and  Az,  which  merge  at  about  110  km 
altitude,  correspond  to  the  normal  and  disturbed  l)-regions , *  respectively 


FIG.  IV-8  MODELS  OF  THE  DAYTIME  IONOSPHERIC  ELECTRON 

RECOMBINATION  COEFFICIENT  ABOVE  80  km  ALTITUDE 


It  is  assumed  that  the  quiescent  »- 
that  the  disturbed  D-region  ion  is 
the  first  assumption  may  not  always 
1  Obi' ) 


region  ion  above  80  km  is  NO1  and 
in  general  a  mixture  of  NO*  and  0' ; 
be  justified  (Poppoff  and  Whitten, 
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This  Model  is  based  on  the  values  of  for  the  various  ionic  species 
given  above,  on  the  Middle  latitude  temperature  profiles  corresponding  to 
quiescent  and  disturbed  conditions  found  by  Spencer,  Brace,  and  Carignan 
(1962),  and  on  the  relative  concentrations  of  ionic  species  as  reported 
by  Taylor  and  Brinton  (1961).  It  was  constructed  under  the  assumption 
that  dissociative  recombination  is  dominant  in  this  altitude  interval  We 
justify  this  by  comparing  the  estimated  value  of  t^(N0+)(  10“7  cm3  sec-1 

with  the  estimated  value  of  the  product  of  the  negative  ion  concentration 
tatio  0  02,  sec  Section  IV-A,,  and  the  ion-ion  recombination  coeffi¬ 
cient  (  2  x  10  8  cm3  sec-1,  see  Section  IV-D);  obviously  a  (uo+)  \  q 

and  we  thus  neglect  the  ion-ion  recombination  processes  in  our  model 
lor  this  reason  the  profiles  also  serve  as  models  or  the  effective 
recombination  coefficient.  The  sharp  rise  in  curve  A,  between  80  and 
85  km  is  a  result  of  our  assumption  that  N0+  is  the  dominant  ion  below 
85  km  under  quiescent  conditions  (Kicolet  and  Aikin,  1960).  If  o*  should 
prove  to  be  dominant  in  this  interval,  curve  A2  will  coincide  with  curve 
A,  below  100  km  It  is  apparent  from  the  dotted  curve  that  for 
k  i  j  _  10  i m  sec  ',  reaction  (It  will  not  have  an  important  bearing 

on  Q()  above  80  km.  The  difference  between  the  curves  corresponding  to 
disturbed  and  quiescent  conditions  at  higher  altitudes  is  due  to  higher 
electron  temperatures  expected  under  disturbed  conditions  (Spencer, 

—  —  ’  Mitia  s  ,1959;  noontime  (quiescent;  model  for  middle 

latitudes  is  also  presented  It  is  evident  that  above  90  km  the  two 
quiescent  models  agree  to  within  a  factor  of  two  The  poor  agreement 
at  lower  altitudes  is  due  in  part  to  (in  our  opinion  Mura's  probable 
over-estimation  of  the  importance  of  electron  attachment  ami  ion-ion 
rccombi nat ion . 

In  the  altitude  range  8(1  to  (.0  km  the  pressur.  rises  to  >  0  1  mm  Hg 
and  the  dissociative  recombination  coefficient s  ol  ionized  molecular 
oxygen  is  expected  to  rise  to  IIP*  cm3  see  (Kasner,  et  al  ,  1961  as  the 

°3  an<l  po,'hapri  °4  U’ns  I’t’i'ome  significant  and  then  dominate  Whether 
or  not  such  clustering  occurs  m  the  nitric  oxide  ion  is  unknown  at 
present  Its  discovery  should  occasion  no  surprise. 
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D. 


Ion-Ion  Recombination 


lo„-lo„  reco.dbl nation  Is  „f  potential  importance  In  the  D-reg,o„ 

Ut,der  b°*h  d"y  *"d  n‘*ht  =™««o„s;  the  altitudes  „  „„,ch  o„.  ™„ 
basin  constdar  1,  arc  or  the  order  ot  70  to  75  to  dun*  ,hc  day  and 
80  to  85  to  at  night .  For  Ionospheric  disturbance.  Witch  ra.ul,  tn  large 
Ionization  rate,  at  low  altitudes  (~  30  to  60  to)  such  as  polar  cap 
event,  and  nuclear  bursts,  ion-ton  reco.bln.tipn  „iu  bo 

recombination  mechanism  and  must  be  Included  in  the  kinetic,  or  the 
ionic  reactions. 

The  ion-ion  neutralization  reactions  are  basically  of  two  types, 
two-body  and  three-bodv.  Typical  examples  are 


Two  Body 


Nj  +  o“ 


02  +  Oo 


n2  +  02 
°2  +  02 


N0+  4  02  .  No  +  o2 


02  -f  0? 


°  +  0  4  o2  (ion-ion  dissociative 
recombination) 


Three  Body 


N2  +  02  4  M  N2  4  02  4  M 


(1) 

(3) 

(3) 

(•1) 

(5) 


where  M  is  an  unspecified  third  body.  Neither  mode  has  been  extensively 
i nvesti Rated. 


Yeung  1058)  measured  the  rate  constants  for  the  process 

^2  ^  1 2  *  1 5  4  I , 


Br - 


*  Br, 


Hr 


(6) 

(7) 
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and  found  values  of  1.47  x  10-8  and  1.85  x  10“8  cm3  sec-1 ,  respectively, 
at  temperatures  of  the  order  of  300 °K .  His  attempt  to  obtain  the  temper¬ 
ature  dependence  of  the  reaction  was  unsuccessful  although  he  did  find 
that  the  rate  constant  decreased  with  increasing  temperature.  In  spite 
of  great  care  in  carrying  out  his  investigations  Yeung  may  have  used  a 
defective  amplifier  in  his  experimental  apparatus  (Sayers,  1962).  Accord¬ 
ing  to  Sayers  the  necessary  correction  under  these  conditions  would  raise 
the  experimental  values  of  the  two  coefficients  to  ~  10~7  cm3  sec-1. 
Unfortunately,  no  experimental  work  has  yet  been  done  on  the  very  impor¬ 
tant  reactions  (1)  to  (4).  Nawrocki  and  Dalgarno  (Nawrocki ,  1961a)  have 
suggested  the  adoption  of  a  mean  atmospheric  ion-ion  recombination  coeffi¬ 
cient  of  a.  =  10-8  cm3  sec-1. 

l 

Many  years  ago  Thomson  obtained  a  semi -empirical  equation  for  the 
three-body  ion-ion  recombination  process 

QT  =  77dD  +  (8) 

where  d^  is  the  effective  »adius  of  coulombic  attraction,  e  is  a  proba¬ 
bility  function,  and  C  is  the  mean  ion  velocity.  More  recently  Massey 
(1952)  showed  that  for  air  at  pressures  <  300  mm  of  Hg  equation  (8)  can 
be  expressed  in  the  more  useful  form 

aT  =  8  x  10" 3  d  T-5/2  (9) 

where  p  is  press  ire  in  mm  of  Hg  and  T  is  temperature  in  °K .  Because  the 
rate  of  two-body  ion-ion  recombination  exceeds  that  of  the  three-body 
mode  at  altitudes  above  45  km,  the  latter  mechanism  is  of  no  importance 
in  the  normal  ionosphere.  It  may,  however,  play  an  important  role  at 
relatively  low  altitudes  during  a  polar  cap  absorption  event  or  during 
a  nuclear  explosion-induced  ionospheric  disturbance. 
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E. 


Diffusion  to  Dust 


Several  years  ago  Bourdeau,  Whipple,  and  Clark  (1959)  measured 

atmospheric  conductivities  f'hnth  n 

(both  P°sl“ve  and  negative)  in  the  altitude 

ranee  35  «  k„  by  of  a  Gcrdien  condenser  by  .  vlkl„E 

J  e,.  Bo.o,  50  k.  the  resu„s  contained  „o  surprises,  b„t  tta  beh,vlo 

01  the  measured  conductivities  was  ni,|.„  , 

33  qUlte  ano^lous  between  50  and  80  km. 
The  positive  ion  conductivity  was  fn,ma  *  , 

nvity  was  found  to  increase  more  slowly  than 
one  -t.it,  norma,, y  _  J 

~  *"  “  ■»  al tl tude  o,  about  30  k. 

I  —  -  co-workers  were  able  «o  tbelr 

.  ’  —or,  ,lth  a  number  densitv 

particle  cm  and  a  particle  radius  of  ^  10-*  Cm  to  b 

.n  ,15  ,e,:i';":;,p::::i ,n 

o!  I”  PrOK"‘0S  "■*  bC“  -  1  nvest lfiators  o,  ,he  ...pro 

»oco,„  ,ZZ  however)  £*£, 

concen t ra t i on  of  dust  »,  ab„u,  „  t,  „  sever.,  orders 
sma  1  lev  than  one  finds  nocm„ce,„  elouds. 

On.  e,„  th,  .ftoc,s  al, „„„  ion  ^ 

T  by  ‘n,'od"«i™  »'  appropriate  diffusion  ten.,, 

the  forr.i 


~a2  C  N„  n 


-here  a  is  the  particle  radius,  No  the  particle  cc 

(t)  i0"  C°nCCntrati0n-  and  C4  the  (i)  ion  rms  velocity.  The  rat^ 
equations  are  written 


d.V 


dt  q  '  V+No  "  V-\  -  a2  CvN„N4 


a2  C+V-  (1) 


d.V 


dt  ~  M  N'e  "  \N’  N  ~  (*  .  ,  +  p  )jf 

1  +  -  'coll  photo'1  - 


a2  N0C  N  -  a2  C  \  \’d 


(*) 


K8 


dN 


dt  "  Q  +  (Ccoll  +  Pphoto)N-  +  6N-  "  QDN+Ne  •  MNe  -  tra2  -  Fa2  C  N0N 


e  u  e 


dN 


dt~  =  lra*  C+NoNJ.  ‘  7Ta‘  C  NdN  -  2-ra2  c  \d  -  t- ,  *  ~  »  -d 

r  +  e  +  e  d  + 


(3) 

77a“  C_  N_  N“  (4) 


dN 


c  -e”°*'e  “  a*  Cd  N+N-  "  6N_  “  C+N+Nd  (5) 


=  Tra*  C_  N0N^  +  rra2  c  x0N  -  2  -2  ~  ■"d-d  ~--d 

assuming  that  every  ion-particle  collision  results  in  charge  capture 
by  the  dust  particle.  The  symbols  have  the  significance 


Ne  =  electron  concentration 

q  =  electron  production  rate 

aD  ”  dissociative  recombination  coefficient 

ai  ~  two-body  ion-ion  recombination  coefficient 

"  specific  electron  attachment  rate 

Pcoll  =  specific  electron  detachment  rate  (col  1  i -m  onal  - 
type)  ' 

photo  '  specific  electron  photo  detachment  rate 

rate  of  electron  detachment  from  negatively 
charged  dust  particles 


Ce  -  electron  rms  velocity 

vd 

concentration  of  positive  (negative)  charges 

residing  on  dust  particles  (electronic  charges 
per  cm3)  K 

Cd  =  dust  Particle  rms  velocity. 

The  third  terms  appearing  on  the  right-hand  sides  of  (4;  and  (5)  can  be 
neglected  because  of  the  very  small  velocity  of  the  dust  particles. 

Inder  equilibrium  conditions  the  concentration  of  positive  charges 
residing  on  dust  particles  is  given  approximately  by 


C  Nn  N 


C  N  +  C  N 
e  e  - 


(6) 
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Assuming  a  temperature  of  20n°K  «  j,  , 

or  200  K,  a  dust  particle  concentration  and  mean 

radius  of  N0  ~  l  cm-3  and  a  ~~  in~s 

10  cm,  respectively,  and  an  electron  con- 

d  °n  °f  Ne  10  Cm  *'  WG  find  that  according  to  equation  (6) 

N+  ^  200  cm'3.  Actually  lt  ls  not  iikely  that  a  dust  parUc]e  of  thig 
size  could  accumulate  more  than  a  few  electronic  charges  because  of 
coulomb  repulsion  of  positive  ions  by  the  positively  charged  dust  particles. 

The  electron  concentration  at  any  altitude  under  conditions  of 
equilibrium  is  obtained  by  adding  equations  (2)  and  (3) 

q  +  ^  =  (l+X)(Vn^  *  V>e  +  ,a3(Ce+Xd  CjM#Me  *  TTa*  C_  N_  N+d 


=  Ne  '  ^  °rdCr  l°  R3in  SOme  lden  of  the  effects  of  diffusion 

to  dust  we  compute  N  at  several 

e  01  altitudes  under  the  assumption  that 

N°  =  1  C'n'1'  Nd  =  1  cm'3>  %  =  lO"7  cm3  sec” 1  ,  and 

=  1°  cm3  sec-,  we  further  assume  that  the  negative  ion  terms  in 
(  ro  *opo„do„,  tho  c„n„nlr„lon> 

The  results  are  shown  in  Table  1V-7. 

It  is  evident  that  particulate  matter  can,  if  present  In  even 
mi nuscu le  quantit.es,  be  very  effective  electron  removal.  view  of 

the  results  of  Volz  and  Goody  (1962)  it  is  doubtful  that  the  mechanism  is 
normal  Is  of  importance  in  the  D-region. 

It  is  easily  shown  that  as  the  electron  concentration  ,s  subs.an- 

tU,Hy  raiSt’°’  PtMhaPS  S‘'V,',al  -gnitude  in  the  case  of  a  high 

antt.de  nuclear  burst,  diffusion  to  dust  becomes  an  even  less  important 

electron  removal  mechanism  unless  thc  dus,  concentration  is  simultaneously 
■ne eased  by  many  orders  of  magnitude.  The  increase  of  the  dust  coucen- 

;;:,0n  l°  N°  10  Cm"S  ‘n  3  l°-km  tUbl  the  proper  dispersal  of 

kg  of  material.  Hence  it  is  not  considered  liKeiy  that  after  such  a 

but  s  l ,  ci  i  f  f u s  i on  of  ol or  l  miK  u.  •»  1  i  v 

°ns  *U1  >>''  important  outside  thc  locale  of 

thc  explosion.  Inside  this  region  it  miv  hr  nr 

M  11  bt  °r  some  significance,  dc - 

pr~.U„B  „po„  ot  dlspr,.,n,  p„.„ru|alc  „iu,r 
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Tabic  IV-7 


ELECTRON  CONCENTRATIONS  WITH  AND  WITHOUT  DIFFUSION  TO  DUST 

(daytime ) 


A 1 1  i  tilde 
(km) 

} 

q 

(cm- 3  sec  1 ) 

*o 

(cm” 3) 

■V 

(wi thout 
di  f fusion ) 
(cm-3) 

Nt. 

(#i th  diffusion) 
(cm"3) 

90 

10-2 

10 

1 

104 

5  x  103 

80 

10'1 

0.1 

i 

10s 

100 

70 

2 

•1  N  10" 3 

i 

103 

1 

60 

30 

10-3 

i 

10 

0.3 

90 

lo-2 

10 

0.01 

10  4 

104 

80 

io-' 

0.1 

0.01 

10 

500 

70 

2 

1  x  10"3 

0.01 

102 

25 

60 

30 

10-3 

0.01 

10 

1 

91 


F '  Effects  of  Chemical  Reactions 
1  •  rnt  roduction 

..  T"°  Ch°rUCaI  comP°s^  ti on  of  the  upper  atmosphere  determines  in  part 

"  t,,0,‘™IC  SPOC‘"  radiation,  „„ 

°f  lh"  »y  ~a„a  of  chame  transfer  and 

It  lmC'Ch"“"  «  -  ionic  sped  os 

~  '"**"*  ™~»*.  «"«ica,  composition  can  i, 

.  *  *“l  •»*«■»  -n  ionospheric  relaxation. 

The  normal  composition  of  the  upper  atmosphere  1,  dependent  not 
I  ;  open  the  Chemical  reaction,  include  photo-reaet. on,  b,„  also  upon 
"*•  ■“  —  For  example,  at 

Tr  “  110  'hC  •'  oaVdcn  atom,  „  0, 

“VCra‘  nocomhinfnp  a,  these  altitudes  after 

formation  hy  phot.-di  „oc,  „„  ,  the,  diffuse  do.„„rd  and  recomb., 

below  100  km  altitude  (Nicolet  1959 

'  1959,  ,  Meanwhile,  oxyBen  molecul. 

U  0  UP'ka'd’  th°  n°l  cffcct  Pci  np  the  establishment  of  diffusi, 
odu...,,rt„„  a,  an  altitude  probably  not  Brc»,cr  than  12„  k., 

UT-  °  '  bC‘°''  '  90  k’’  «-  -Chine  a,  „„  they 

°  '  or  diffuse  ,„,o  the  region  and  „  have  a  condition 

t,'  77“'  ThC  -  *»'ch  turbulent  mlx.nff 

e,„d, b-"“  * . .  'no.  Photo-che.ical 

.  a  ~  lh"  "a“""C  «  *>*  -»»•*»  opuilibrium 

‘  Inactions  o,  . . .  o,  dissociation  and  scale  belch, 

;;; ; . .  fi*  •-  >»  ™  - .* 

; ;;  ™i'ibr>m -  — *•* «... « 

cqui  1 1  hi  itiif  (Nicoict,  19(50 


ine 

les 
t  ve 


II  -do  system  is  be  maintained  photo-che.ioal  eunilihrlum 

. . . 

/  ‘  1 11  '*  . .  2  to  3  day.;  „„s  rcoulre^n,  is  due  to 

r“°  “'VBC0  ’ hroughout  „,e  day 
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Recently  several  model  atmospheres  have  been  developed  on  the 
bases  of  satellite  drag  measurements,  measurements  of  the  absorption  of 
solar  radiation,  and  estimates  of  the  heating  effect  The  one  presented 
in  Figure  IV-9  is  due  to  Norton,  Van  Zandt ,  and  Denison  (1962)  for  alti¬ 
tudes  above  100  km  and  the  Handbook  of  Geophysics  I960  for  altitudes 
below  100  km,  whil.  t  ha  t  pies.ntcd  in  Figure  IV-10  is  taken  from  a  paper 
b\  Barth  1901  and  is  based  on  the  work  of  Ko 1 lman-Bi jl  (1961  and 
Nicolet  1959  The  foimer  Figure  IV-9)  is  probably  the  more  reliable 
since  it  is  based  in  part  on  measurements  of  the  absorption  of  solar 
ladiation  The  Handbook  of  Geophysics  model  atmosphere  was  selected 
foi  the  lower  altitudes  in  Figure  IV-9  because  it  joins  smoothly  with 


TEMPERATURE  —  *  K 

**  i*?4  2t 

FIG.  IV-9  ATMOSPHERIC  NEUTRAL  CONSTITUENTS  AND  TEMPERATURE, 

60  TO  180  km  ALTITUDE  Norton,  et  ol,  1962  ond  Hondbook  of  Geophysics,  1960) 
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*'  »•!«  it 


FIG.  IV- TO 


ATMOSPHERIC  NEUTRAL  CONSTITUENTS  AND  T 
60  TO  160  km  ALTITUDE  (Borth,  1961) 


EMPERATURE. 


the  model  atmosphere  of  Norton,  ct  al  at  100  km, 
lo  which  the  latter  was  extended 


the  lowest  altitude 


Basic  Chemical  Reaction* 


Barth  1901J  has  thoroughly  discussed  the  chemic 
the  upper  atmosphere  which  control  the  concentrations 
molecular  oxygen  and  nitrogen,  nitric  oxide,  nitrogen 
07000 '  Thest*  tactions  can  be  summarized  as  follows 


al  reactions  in 
of  atomic  and 
dioxide,  and 


N  +  N  +  M  N2  +  M 

N  +  0  +  M  -  NO  +  M 

0+0+  M-  02+M 

N  +  NO  N2  +  0 

+  02  -+  NO  +  0 
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N 


o  +  o2+m-»  o3+m 

(6) 

o  +  03  -»  02  +  02 

(7) 

NO  +  0  +  M  -4  N02  +  M 

(8) 

N02  +  0  -  NO  +  02 

(9) 

With  the  corresponding  rate  constants  obtained  by  various  investigators 
shown  in  Table  I V-8  In  addition,  Table  IV-8  contains  in  the  right-hand 


RATK  CONSTANTS  CORRESPONDING  TO  PROCESSES  (1  TO  (9 


Ra  te 

Constant 


k, 


kc 


Barth  (1961  / 


7  x  10" 3  3  err6  sec-1  (300°K 


1.5  x  10“ 32  cm'  sec-1  (300°K) 
«  0  x  10" 3 3  cmc  sec"  1  (300 °K 


«  3  x  10' 11  cm3  sec"1  300°K 

_  3 1  oo 

3  3  x  10” 1 2  e  T  cm3  sec-1 


6  8  x  10- 


300 

T 


.3000 

T 


cmG  sec-1 


5.0  \  10” 1 1 
5.3  s  l()“3J  cm6  sec 


3  5  x  10” 1 2  cm3  sec” * 


cm3  sec-1 


Hurt  nor  and  Bnulkniglit 


1 

.7 

x  lo¬ 

30  j 

T-i 

T-* 

T-i 

-4  6 

cm”  sec 

1  . 

(M 

N2). 

2 

X 

ur31 

cm 45 

see” 1 ; 

(M 

NJ* 

3 

X 

10"31 

cm6 

sec  “* , 

(M  p 

N2.  NJ 

5 

X 

10-3* 

cm6 

sec”  1  • 

6 

X 

10*  33 

T*i 

T~  * 
T-i 

cm4. 

sec”  1  ; 

(M 

— 

02)t 

2 

X 

lO-3' 

cm6 

see" 1 ; 

(M 

= 

0)t 

1 

X 

10“32 

cmB 

) 

sec-1 ; 

(M 

t 

o 

c 

M 

2.3  x  10- 
1  x  10-' 3 
1  x  10"  34 


300°K 


300°K 


7  x  ID-  1 2 
5  x  10“ 3 2  cm-  sec' 


T  cm3  sec” 1 • 

_  66  oo 

A  RT  , 

T  e  cm3  sec" 
-  LiP 
T 

e  cm1’  sec  -  1  • 

_  I  OOP 

T 

'  cm 3  my’1  ■ 


2  5  x  10“ 1 2  cm 3  sir-1 


tBortncr  and  Baulknight  1961 
‘Bortner  and  Baulknight  1962 

*Foi  specific  references  see  Barth's  v1901  papei 
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column  a  set  of  rate  constants  proposed  by  Bortner  and  Baulknight  (1961, 
1962)  which  are  somewhat  different  from  those  of  Barth. 


3 .  Chemistry  of  the  Normal  Atmosphere 


In  order  to  write  and  solve  differential  equations  for  the  rates 
of  change  of  the  concentrations  of  the  various  species,  one  must  know 
the  rates  of  dissociation  of  nitrogen  and  oxygen.  The  dominant  mechanism 
for  the  formation  of  atomic  oxygen  in  the  upper  atmosphere  is  photo¬ 
dissociation  of  02  by  solar  radiation  in  the  Schumann-Runge  continuum 
(1300-1750  A).  The  photo-dissociation  cross  section  in  this  region  is 
of  the  order  of  5  x  10  **  cm  (Watanabe,  1958)  and  the  solar  photon 
flux  ~  5  x  1013  cm2  sec  1  (Mitra,  1952).  Using  the  concentration  of  02 
given  in  Figure  IV-10,  and  neglecting  attenuation  of  the  photon  flux  by 
absorption,  one  finds  an  approximate  value  of  the  rate  of  dissociation 
of  Oj  of  5  x  10*  cm'3  sec'1. 


On  the  other  hand  photo-dissociation  of  N,  is  a  very  weak  process 
The  N2  molecule  is  apparently  dissociated  by  absorption  in  the  Lyman- 
Bi rge-Hopf ield  band  at  1226  A,  for  which  the  estimated  cross  section 
is  about  10  21  cm2  Since  the  photon  flux  is  ~  109  cm-2  sec'1  (Barth, 

1961)  we  obtain  a  photo-dissociation  rate  at  100  km  of  only  10  cm-3  sec-1. 
A  process  which  promises  a  much  larger  rate  of  formation  of  nitrogen 
atoms  is  the  dissociative  recombination  of  N2  or  N0+ .  In  order  to 
obtain  an  estimate  of  the  effective  dissociation  rate  of  NO*  at  100  km 
we  estimate  that  [e]  =  2  x  10s  cm-3, (NO*]  *  10s  cm'3  and  ^(NO*)  5  x  10'* 

cm3  sec-1.  Then  the  rate  of  formation  of  nitrogen  atoms  from  the  process 


NO1  »  e  .  N  +  0  10 ) 

‘  *j>  N0  f'°  i‘‘J  1()3  3  sec-1,  a  much  larger  value  than  that 

arising  from  photo-dlssoclatlon  of  N2 


It  is  in  principle  now  possible  to  compute  the  chemical  composition 
of  the  upper  atmosphere  if  one  knows  how  to  handle  the  diffusion  problem 
In  practice  this  is  a  formidable  problem  and  Barth  contented  himself 
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with  starting  with  the  initial  (daytime)  conditions  shown  in  Figure  IV-10 
and  computing  the  concentrations  of  various  constituents  at  various  times 
after  sunset.  It  was  found  that  the  concentration  of  atomic  oxygen  was 
quite  constant  during  the  night  for  altitudes  above  70  km  while  the 
concentration  of  atomic  nitrogen  remained  constant  at  altitudes  of 
80-120  km  but  decreased  rapidly  during  the  night  at  altitudes  above  or 
below  this  range.  The  decrease  of  Oj  and  [N]  during  the  period  between 
sunset  and  one  hour  thereafter  are  shown  in  Figure  IV-11,  under  the 


FIG.  IV-U  CHANGES  IN  CONCENTRATION  OF  0,  N,  AND  NO  AFTER  SUNSET 
AT  ALTITUDES  OT  60,  80,  AND  100  km 


assumption  that  during  the  daytime  Nj  =  10s  cm” 3 

Formation  and  destruction  of  nitric  oxide  result  primarily  from 
reactions  (5)  and  (4)  respectively.  These  processes  togethe'r  with 
diffusion  yield  the  daytime  NO  concentration  curves  shown  in  Figure  IV-10. 
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er  sunset  is  show  in 


The  decrease  of  [ NO J  during  the  first  hour  aft 
Figure  IV-11. 


The  composition  of  the  upper  atmosphere  has,  of  course,  a  pronounced 
effect  on  ionospheric  structure  since  it  determines  to  a  large  extent 
the  nature  of  the  ions  present,  the  ionization  rate,  and  the  rate  of 
electron  removal.  However,  this  composition  is  nearly  constant  under 
normal  conditions,  the  concentrations  of  a  few  of  the  minor  constituents 
at  certain  altitudes  decreasing  after  sunset.  The  only  change  of  this 
kind  which  may  significantly  alter  the  electron  removal  rates  is  the 
combination  of  atomic  and  molecular  oxygen  to  form  ozone  at  altitudes 
below  70  km.  Specifically,  such  a  reduction  in  the  atomic  oxygen 
concentration  will  result  in  a  smaller  rate  of  associative  detachment 
of  electrons  from  negative  ions  (e.g.,  oj  at  night  than  in  the  daytime. 

The  change  of  the  solar  spectrum  between  maximum  and  minimum  of 
•solar  activity  undoubtedly  affects  the  composition  of  the  upper  atmosphere 
but  the  details  arc  not  at  all  well  understood. 

'*  •  Effects  of  lij  st  u rbanccs 

Abnormal  solar  radiation  of  short  duration  which  is  capable  of 
producing  changes  in  the  composition  of  the  upper  atmosphere  may  be 
broadly  classified  as  electromagnetic  (X-rays  and  ultraviolet  radiation) 
and  Particle  ,sola,  cosmic  rays  The  increased  rate  or  dissociation 
<>f  0,  arising  from  tin  enhanced  radiation  is  certainly  not  greater  than 
103  cm’3  sec’1  ..s  opposed  to  a  normal  rate  of  dissociation  of  p2  at 
90  km  of  10“  cm  3  sec-1,  the  enchancement  is  therefore  negligible. 

On  the  other  hand  the  formation  of  nitrogen  atoms  proceeds  via  the 
processes 

•V  3  "  •  N  +  N  (11) 

NO'  +  e 


N  i  0 


10 


rather  than  by  direct  dissociation  of  N2 .  It  is  demonstrated  in  the 
section  on  ion-atom  interchange  that  charge  transfer  to  02  by  N2+  is 

rapid  a  pt  ’cess  that  N2  is  a  negligible  ionic  component  of  the  lower 


9X 


I 


ionosphere,  although  it  is  produced  at  a  faster  rate  than  any  other 
ion.  Undoubtedly  some  N0+  is  formed  by  the  ion-atom  interchange  process 

N2  +  0  — »  NO  +  N 


which  is  rather  slow  D.  R.  Bates,  private  communication)  as  compared 
to  other  processes  of  similar  type.  Assuming  that  [no*  u  as  larRe 
as  0.1  02  ,  we  obtain  a  production  rate  of  N  atoms  of  -  10  cm"3  sec'1 
and  an  enhancement  in  the  N  concentration  of  <  104  cm  3  which  is  less 
than  the  normal  concentration  thought  to  arise  via  diffusion  from  higher 
altitudes.  Even  if  the  change  in  Nj  were  significant,  the  effect  on 
ionospheric  relaxation  would  still  he  negligible. 

In  the  case  of  a  high  altitude  nuclear  burst,  the  molecular  species 
are  almost  completely  dissociated  within  the  fireball.  Just  outside 
the  fireball  region  ultraviolet  radiation  in  the  range  1300-1700  A  is 
sufficiently  intense  to  produce  virtually  complete  dissociation. 

However,  at  altitudes  below  90  km  the  attenuation  of  this  radiation  by 
absorption  is  so  large  that  it  has  a  range  of  .  1  km,  assuming  an 
absorption  cross  section  of  10—  cm*  (Watanabe,  1958).  At  an  altitude 
of  100  km,  dissociation  of  02  is  already  so  great  that  nearly  complete 
dissociation  of  the  residual  02  will  no,  greatly  affect  the  concentration 
of  atomic  oxygen  although  the  concentration  of  02  may  be  decreased  by 
ns  much  as  an  order  of  magnitude  until  such  time  that  recombination  and 
diffusion  can  restore  the  system  to  equilibrium. 


If  we  assume  that 
is  of  the  order  of  ]()8 
concentration  of  atomic 
10'°  cm-3  which  may  in 
an  order  of  magnitude; 
intensity  measurements 
not  sufficiently  great  , 
loss  mechanism. 


lh.  mean  electron  density  at  an  altitude  of  loo  km 
cm  3  over  a  period  ol  at  least  one  hour,  the 
nitrogen  may  be  raised  to  a  value  as  large  as 
turn  increase  the  concentration  of  NO  by  perhaps 
the  increase  in  N  may  be  detectable  by  nirglow 
VN°  bands  These  changes  in  composition  are 
however ,  to  significantly  affect  the  electron 
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